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ABSTRACT 

We present a study of the mechanical power generated by both winds and jets across the black 
hole mass scale. We begin with the study of ionized X-ray winds and present a uniform analysis 
using Chandra grating spectra. The high quality grating spectra facilitate the characterization of the 
outflow velocity, ionization and column density of the absorbing gas. We find that the kinetic power 
of the winds, derived from these observed quantities, scales with increasing bolometric luminosity as 
log(L w i n( j : 42/C„) = (1.58 ± 0.07) log(LBoi,42) — (3.19 ± 0.19). This suggests that supermassive black 
holes may be more efficient than stellar-mass black holes in launching winds, per unit filling factor, 
C' v . If the BHB and AGN samples are fit individually, the slopes flatten to a BHB = 0.91 ± 0.31 
and a AGN = 0.63 ± 0.30 (formally consistent within errors). The broad fit and individual fits both 
characterize the data fairly well, and the possibility of common slopes may point to common driving 
mechanisms across the mass scale. For comparison, we examine jet production, estimating jet power 
based on the energy required to inflate local bubbles. The jet relation is log(Lj c t,42) = (1-18 ± 
0.24) log(LBondi,42) — (0.96 ± 0.43). The energetics of the bubble associated with Cygnus X-l are 
particularly difficult to determine, and the bubble could be a background supernova remnant. If we 
exclude Cygnus X-l from our fits, then the jets follow a relation consistent with the winds, but with a 
higher intercept, log(Lj e t,42) = (1-34 ± 0.50) log(LB on di,42) — (0.80 ± 0.82). The formal consistency in 
the wind and jet scaling relations, when assuming Lb i and LBondi are both proxies for mass accretion 
rate, suggests that a common launching mechanism may drive both flows; magnetic processes, such 
as magneto-hydrodynamics and magnetocentrifugal forces, are viable possibilities. We also examine 
winds that are moving at especially high velocities, v > 0.01c. These ultra-fast outflows tend to 
resemble the jets more than the winds in terms of outflow power, indicating we may be observing 
a regime in which winds become jets. A transition at approximately Lb i ~ 10~ 2 LEdd is apparent 
when outflow power is plotted versus Eddington fraction. At low Eddington fractions, the jet power 
is dominant, and at high Eddington fractions, the wind power is dominant. This study allows for the 
total power from black hole accretion, both mechanical and radiative, to be characterized in a simple 
manner and suggests possible connections between winds and jets. X-ray wind data and jet cavity 
data will enable stronger tests. 



1. INTRODUCTION 

Both winds and jets are thought to be driven by ac- 
cretion disks; jets may be launched from the innermost 
regions, while wi nds may originate further out in the ac - 
cretion disk (e.g. JBlandford fe Pavndll^82l:lProgal[2003l ). 
It remains to be seen just how winds and jets relate. Are 
winds and jets driven by similar mechanisms? Does one 
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quench the other? What role does the geometry of the 
magnetic field lines play? One might expect the same 
physical launching mechanisms across the mass scale. 
Jet production is predominantly ascribed to magneto- 
hydrodynamics (MHD) across all mass sc ales; whether it 
is purely through MHD i n the disk (e.g.. lLovelacdll976T ; 
Blandford fe Pavnelll982f), or through t he disk and black 
hole (e.g., iBlandford fe Znaiekl [19771 : iKrolik fc Hawlevl 
I2010T) . However, it has yet to be shown how wind prop- 
erties scale with mass. Does radiation driving, thermal 
driving or magnetic processes drive these winds? 

"Warm-absorbing" winds detected as X-ray absorp- 
tion features are seen in up to 50 % of Active 

Gala ctic Nuclei (AGN) (jRevnoldsl 119971 : iGeorge et all 
119981 ) and in the soft sp ectral state in stellar-mass 
black holes dMiller et aft l2006alE 120081: lUeda et all 
2009| :lNeilsen fc Ledl2009t King et al]l2012bt iPo^ti et alJ 



2012f h These blue-shifted absorption features are highly 
ionized and can potentially probe regions close to the 
black hole and basic disk physics. 

The winds observed in both black hole binaries (BHB) 
and supermassive black holes (SMBH) are capable of 
removing enormous amounts of ma terial, even exceed- 
ing the mass acc retion rates (e.g., iBlustin et al.1 120051 : 
iKing et aill2012bl ). They are typically wide angle out- 
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flows moving at a few hundreds of km s _1 . Jets, on 
the other hand, are highly collimated and have a much 
higher outflowing velocity, i.e. near the speed of light. 
Consequently, the mechanical energy in these jets can 
be much higher than in winds. In addition, using only 
the radiative luminosity in jets severely underestimates 
the power released in these systems because a major- 
ity of the energy is mechanical (e.g.. iGallo et al.l [20051: 
Churazov et alj|2005t lAllen et ali 2006; Merloni fe Hemzl 



20071 ). More importantly, strong mechanical feedback 
from a black hole can have a significant impact on 
its surroundings, includin g galactic formation, struc- 
ture and co-evolution (e.g..lHaehnelt fe Kauffmannll2000l: 
Churazov et al.ll2002HCroton et all 120061: lOstriker et al.1 



2010: iGaspari et al. [ 120111: lFabianll2012D . 



On the low end of the black hole mass scale, the driv- 
ing mechanisms of BHB X-ray winds are generally as- 
cribed to either thermal pressure or magnetic mecha- 
nisms (|Begelman et alJ[T98l iWoods et al1[l99l |P~" 



rog, 



a 



|2003|) . Absorption features of highly ionized Fe XXV and 
Fe XXVI are the most commonly detected, implying the 
ionization parameters of the gas in BHB are ver y high , 
i.e. log£ > 3 ([Miller et all I2006bl: King et al.l l2012bl 



iPonti et al.l 12012 ). At such high ionization states, line 
driving from the radiation field is inef ficient at accelerat- 
ing t he winds to high velocities (e.g., iProga fe Kallmanl 
2002). On the high mass end, X-ray winds from SMBH 
systems span a larger range in ionization (0 < log£ < 5). 
Therefore, not only are Compton heating and magnetic 
mechanisms plausible driving mechanisms, but radiation 
pressure at low ionizations is also a plausible driving 
mechanism (jProga et al.ll2000h . 

In this paper, we begin to examine the mechanical out- 
flow power released in both winds and jets. In addition, 
we examine the outflow power across the mass scale, in- 
cluding stellar-mass and supermassive black holes. This 
will permit a complete characterization of the output of 
black holes: both radiative and mechanical. Character- 
ization of the mechanical power can be particularly im- 
portant on larger scales with respect to AGN feedback. 
AGN outflows, both jets and winds, may be responsi- 
ble for shaping their environment, whether providing a 
source of hot ionized gas or by influencing the stellar ve- 
locity dispersion as ev idenced by the M — a relation (e.g., 
IGultekin et al.l [2009bf l . These outflows may also play a 
vital role in the g rowth of black holes (e.g., Hopkin s et al.l 
[2^?j5llLoeD[2005r ). 

We describe how we derive the samples in Section [2l 
Then, we estimate the kinetic power generated by winds 
and jets in these systems in Sectional Next, in Section 
IH we describe the results for both the wind and jet re- 
lations, while in Section [5] we present the conclusions as 
well as context for this study. We assumed Hq — 70 km 
s _1 Mpc -1 , = 0.7, Qm = 0.3 throughout this work. 
All errors are la statistical uncertainties unless otherwise 
stated. 

2. THE MECHANICAL OUTFLOW SAMPLE 

2.1. X-ray Winds 

The central goal in assembling this wind sample is to 
create a uniform, unbiased and cohesive set of standards 
to ensure high quality spectra and rigorous results. These 
requirements are: (1) blue-shifted X-ray absorption fea- 



tures, (2) Chandra grating spectra, (3) photoionization 
modeling, (4) at least a 3cr significant detection, and (5) 
velocity outflows of less than 3,000 km s . 

Wind and jet launching mechanisms are of prime inter- 
est in this analysis, so it is important to probe the wind 
regions that are likely to be closest to the base of the 
flow, and closest to the black hole. This points to a com- 
parison of X-ray winds across the mass scale. Pragmatic 
considerations also make this the only consistent com- 
parison that can be drawn. In AGN, UV radiation pres- 
sure may help to accelerate winds (but radiation pressure 
may not be sufficient to lift gas off of the disk in the first 
place); however, UV winds are not observed in BHB's, 
likely owing to the high ionization parameters found in 
these winds (3 < log£ < 5). Fortuitously, restricting our 
analysis to X-ray winds in AGN captures the bulk of the 
mass outflow rate. Prior treatments of "warm absorbers" 
in AGN have found that the mass outflow rate scales 
with the ioniz ation parameter estimated in different com- 
pone nts (e.g.. iBlustin et all 120051 : iCrenshaw fe Kraemerl 
l2012ll . We also require that observations of the X-ray 
winds be made with the gratings spectrometers aboard 
Chandra. Although the XMM-Newton reflection grat- 
ing spectrometer (RGS) is similar in many respects, it 
covers a narrower energy range, and it has a lower spec- 
tral resolution. The higher resolution of the Chandra 
gratings means that it is more sensitive to lines that are 
intrinsically narrow and weak, because more line flux is 
concentrated in fewer spectral bins. Thus, this selection 
criterion serves the aim of not biasing our result against 
weak or slow X-ray winds. Similarly, we do not consider 
CCD spectra of X-ray winds in our initial analysis since 
the modest resolution of such data inhibits the detection 
of weak lines and modest velocity shifts. 

In addition, we require that the winds be observed as 
blue-shifted absorption features with respect to the host. 
A significance of at least 3<r for each component after 
fiducial fits with a canonical Galactic absorption and a 
power-law is also required. Further, to quantify these 
particular features, we demand that self-consistent pho- 
toionization modeling be performed to determine the gas 
parameters: the outflowing velocity, v ou t, and the ioniza- 
tion parameter, £. The ionization parameter is defined 
as £ = Li on (nr 2 ) -1 , where Li on is the ionizing luminos- 
ity, n is the density, and r is the distance from the warm 
absorber to the ionizing source. There are also references 
in the literature to the ionization parameter U, which is 
defined as U = Q(47mr 2 c) _1 , where Q is the number 
of ionizing photons and c is the speed of light. We use 
nr 2 = Q(4ttUc)~ 1 — iion£ _1 to convert to the ionization 
parameter £. It should be noted that £ depends on the X- 
ray continuum which is well constrained by observations, 
where as U depends on the number of ionizing photons, 
which is very model dependent. The physical charac- 
terization of ionization in the spectra relies on multiple 
lines to determine the column density, ionization state 
and velocity shifts of the wind components. We include 
fits published in the literature. A majority of those fits 
were obtained with XSTAR grids (|Bautista fe Kallmanl 
|2001|) . which mainly use £ to characterize the ionization 
of the gas. How ever, some fits were also obtained with 
CLOUDY grids (jFerland et al.l |1998f). which mainly use 
U to characterize the ionization of the gas. 
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Finally, we restrict the outflowing velocity to be less 
than 3,000 km s" 1 (0.01c). We reserve outflows with 
a velocity faster than this for a sub-sample of tentative 
ultra- fast outflows. As these fast winds approach such 
velocities they may resemble jets more so than the typical 
lower velocity outflows. 

As a result of these criteria, we select thirteen AGN 
and ten BHB observations. Table Q] lists the ionization, 
velocity and kinetic power for each observed outflowing 
component that is used in this work. Table [5] reports 
the total kinetic luminosity from each observation, sum- 
ming over all the outflowing components. Figure [1] plots 
the observed velocity shift as compared to the ioniza- 
tion parameter, and Figure [2] plots each individual ab- 
sorption component's kinetic luminosity per filling fac- 
tor, C v , versus bolomctric luminosity. We note here that 
our sample may suffer from a few selection biases. The 
first being that we may be biasing ourselves to high lumi- 
nosity sources for which we get the best signal to noise. 
However, our sample does not necessarily include just 
the brightest sources, but also sources that are relatively 
faint and have long exposures in order to increase their 
signal to noise. (See the following sections for details 
of particular observations.) Therefore, a luminosity bias 
should not play a major role in our data. In addition, as 
will be shown in Section [31 the kinetic wind luminosity 
does not directly depend on the column density, i.e. sig- 
nal to noise, but instead on the velocity and ionization 
of the gas. Therefore, as long as a significant detection is 
made, the depth of the absorption features will not serve 
to bias our samples. Another bias could be the exclusion 
of XMM-Newton data because of its poorer resolution. 
However, as noted previously, the lower resolution of the 
RGS would serve to bias the sample against weak fea- 
tures. Further caveats are discussed in Section IS~3l 

2.2. Supermassive Black Holes 

This sample, se en in Table [H is predom inately derived 
from the work by iMcKernan et all ((2007), which is com- 
prised of 15 nearby AGN that were all observed with 
Chandra High Energy Transition Grating Spectrograph 
(HETGS) before July 1, 2003. The summed MEG first- 
order spectra were used. For fu rther details of the data 
red uction, we refer the re ader to Mc Kernan et al.l ([2007D 
and lYaaoob et al.l (|2003ft . 

Owing to the observational selection criteria, only 7 of 
the 15 AGN show statistically significant blue-shifted ab- 
sorption features. IMcKernan et al.l (|2007D perform a uni- 
form analysis with XSTAR models, which characterizes 
the absorption and emission features seen in the sam- 
pled spectra. The XSTAR models were generated as- 
suming an individual SED for each AGN. Each grid of 
models had an assumed density of n e — 10 8 cm -3 and 
a turbulent velocity of 170 km s _1 . These AGN span 
a redshift range of z=0.003 to z = 0.046, a mass range 
of 6xl0 5 — 10 8 M , and a range of environments from 
centers of clusters to field galaxies. 

A majority of these AGN had more than one out- 
flowing component. See Table [TJ These components 
were separated not only in velocity space but in ioniza- 
tion parameter as well. In these cases, we took the sum 
of all the components to evaluate the mass outflow rate, 
M W i n d, and kinetic wind luminosity, L w i n d, so as to com- 



pare the total out-flowing material and consequently the 
total power generated by these winds. See Section [3] for 
details. The values for the kinetic wind power and bolo- 
metric luminosities are given in Table [1] for individual 
components and Table [2] for the summed components. 
We have assumed a covering fraction, SI, of O = 2tt, and 
given the results per filling factor, C v . 

In addition, six other AGN observations are included 
in our sample: NGC 4051, NGC 4593, Mkn 509, IRAS 
18325-5926, NGC 4151 and Mrk 290. These are 
all relatively nearby AGN, z = 0.002 - 0.034, with 
comparable masses spanning M = 2 - 160 xl0 6 M Q . 
King et all ([2012afl report observations of the nearby, 
Seyfert-1 AGN, NGC 4051 which show evidence of warm- 
absorbers. They coadded 12 Chandra HEG and MEG 
spectra from November 2008 for a total e xposu re time 
of 308 ksec. Although, IMcKernan et al.l (|2007h report 
detections o f warm -absorbers in NGC 4051, the work by 
King et all (I2012aft uses a different data set which was 
observed over 5 years later. King et ail ([2012aD use both 
XSTAR and Cloudy photoionization codes to create grids 
of models to fit to the data. By modeling the spectra 
with two separate photoionization codes, they were able 
to determine that three different ionization components 
were required by the data, independent of the model 
used. For this analysis, we use the three components 
that were determined using the XSTAR photoionization 
grids, which span a wide range in ionization and velocity, 
i.e. log£ = 1, 3.4, & 4.5 and v = 400, 630, k 680 

We also used the work by Steenbrugge et al. (2003), 
who observed NGC 4593, both with Chandra LETGS 
and XMM-Newton. We used only the 108 ksec LETGS 
spectrum that was analyzed with an XSTAR model. 
The XMM-Newton observation, taken 7 months later, 
could not constrain the absorption component ioniza- 
tion; however, it was consistent with the LETGS ob- 
servation. The statistically significant absorption com- 
ponent in the LETGS observation had an ionization of 
log£ = 2.61±0.09 with an outflowing velocity of 400±121 
k m s -1 . 

lEbrero et all (|2011h also used the Chandra LETGS to 
observe the AGN Markarian 509. This 180 ksec ob- 
servation was modeled with xabs, created using Cloudy. 
The fit to the spectra resulted in 3 components, one of 
which was significant at the 3cr confidence level. This 
had an ionization of log£ = 2.26±0.07 and a velocity of 
v = 196^ 8 3 k m s _1 , and is included in our sample. 

In addition, we used the work bv lZhang et al.l ([2011a[) . 
who observed IRAS 18325-5926 . The authors coadd 
two exposures from March 2003 to get a total Chandra 
HETG spectrum of 108 ksec. Using a grid of XSTAR 
models, they find two outflowing absorption components 
with typical warm absorbing parameters of log £ = 1.58 ± 
0.09 and log£ = 2.35 ±0.25 and v out = 340±110 km s" 1 
and Vput = 4 60±2 20 km s _1 , respectively. We note that 
iMocz et al.l ([201 ll ) also get results consistent with the 
h igher velocity compone nt . 

iKraemer et al.l ((2005) describe an analysis of NGC 
4151, a nearby Seyfert-2 (it is more like a Seyfert-1 in 
X-rays) , using two coadded Chandra HETG spectra from 
May 2002. The total exposure time was 250 ksec. The fo- 
cus of their work was to describe the absorption features 
seen in both the X-ray and UV using Cloudy models. 
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lKraemer"e t al. ( 20Q5|) do not quote the significance of the 
detection of these features; however, the individual lines 
are detected at several times their minus-side errors, and 
modeled together, should be quite significant. The X-ray 
absorption components taken from lKraemer et"aTI (j2005[ ) 
that are included in our study have ionization parame- 
ters of \ogU ps —0.27 and \ogU ps 1.05 and outflowing 
velocities of v out ~ 500 km s _1 . 

Finally, t he last AGN i n our s ample, Markarian 290, is 
taken from lZhang et all (j2011b[ ). We note these authors 
use both Chandra and XMM-Newton grating spectra in 
their analysis. However, we only include the Chandra 
HETG spectra in our analysis. This particular spectrum 
had three coadded observations giving an exposure time 
of 166 ksec. The two XSTAR grid components that were 
statistically significant spanned a range of ionization 
states from log £ = 1.62 + 0.15 to log £ = 2.45±0.04 and a 
range of velocities from w out = 450±30 to u ou t = 540±150 
km s , respectively. 

2.3. Stellar-Mass Black Holes 

To the greatest extent possible, values in the literature 
were used to estimate the kinetic power in winds ob- 
served in stellar-mass black holes. As with supermassive 
black hole winds, only observations obtained with high 
resolution gratings spectra and some level of photoion- 
ization modeling are included. As a result of these con- 
siderations, only Chandra /HETG spectra were selected. 
In total, ten observations from four stellar-mass black 
holes and black hole candidates are included in our ini- 
tial analysis. It should be noted that many more HETG 
observations are available - GRS 1915+105 has been ob- 
served regularly - but we have only included prominent 
low- and high-flux spectra from GRS 1915+105 in order 
to represent the properties of the group rather than just 
one source. The winds in X-ray binaries have only been 
detected in the "high/soft" X-ray spectral state (e.g., 
IPonti et al.H2012D . We do not include upper limits when 
the sources are in the "low/hard" state. See Section 5.2 
for more discussion of spectral state dependence. 

GRO J1655 40: This is a famous and recur- 
rent transient, and its mass and distance are well- 
determined (M = 7 . ± . 2 M w and d = 3.2 ± 0.2 kpc 
IHiellming fc Rupenl [19951 IQrosz fe BaTlvnl [l997h . The 
presence of density-sensitive Fe XXII absorption lines in 
a Chandra spectrum (ObsID 5461) of GRO J1655-40 
allowed for direct constraints on the density of the dis k 
wind in this source (log(n) = 14: IMiller et "all 120081) . 
Fits with an independent photoionization code, grids of 
Cl oudy models, and gr ids of XSTAR models are reported 
in IMiller et alj ([20081 ). Based on that work, we have 
used values of L = 5.0 ± 0.1 x 10 37 ergs s _1 , an out- 
flow velocity of v = 500 km s _1 , log(£) = 4.9 ± 0.1, and 
O = 2.5 in estimating the kinetic power of the wind in 
GRO J1655-40. The wind observed in GRO J1655-40 is 
particularly complex, and no single velocity characterizes 
all of the lines observed; v = 500 k m s -1 is a represen - 
tative value that is used throughout Miller et al.1 ([20081 ) 
because it captures the outflow well. The resulting ki- 
netic power i s bro adly consistent with numbers given in 
IMiller et al.l (I2008D estimated usin g the wind density. We 
note that iNeilsen fe Homanl (|2012t ) also model a Chandra 
observation (ObsID 5460) that was made a month prior 
to this observation (ObsID 5461). We find evidence for 



an absorption feature at 6.97 keV, which is consistent 
with no outflow, and therefore do not include it in this 
analysis. 

H 1743—322: This source is also a well-known 
and recurrent transient (see, e.g.. IHoman et al.l 120051; 
Corbel et aIJ200alMifler et al.112006 bt IMiller- Jones et al l 
20121) . The high column density along the line of sight 
to H 1743—322 has perm itted the detection of a counter- 
part ([Steeghs et al.ll200"3l) but has not permitted a radial 
velocity measurement. In t his work, we have assumed 
a distance of 8.5+0.8 kpc (Ste iner et al.l [20121) . and a 
fiducial mass of 10 M Q . During its 2003 outburst, a disk 
wind wa s clearly detected in thwo Chandra/UETG obser- 
vations (Mill eret al.l 1200 6b?): parameters obtained from 
photoionization modeling of each spectrum are used in 
this analysis. 

In particular, the broa dband X-ray spectral fits in table 
2 of IMiller et all (I2006bl) were used to derive bolometric 
lum inosities, and the ion ization parameters given in table 
5 m IMiller et al.l ([2006b) are used. The first observation 
included has an outflowing velocity of 670+170 km s _1 , 
while the second observation included has an outflowing 
velocity of 340+170 km s" 1 . The ionization of the two 
observations are roughly the same at log£ = 5.5 ±0.1 
and log£ = 5.6 ± 0.1, respectively. Less is known about 
the binary parameters of H 1743—322 than e.g. GRO 
J1655— 40, and a larger value of ft — 2ir was adopted in 
calculating the photoionization models (the package use d 
was an update of the code described in [Raymond 1993). 

GRS 1915+105: This is an extremely well-known mi- 
croquasar. The mass of the black hole and its distance 
have been determin ed (M = 14 ± 4 M ; d = 12.5 kpc; 
IGreiner et all [2001). A long Chandra/ RET GS observa- 
tion of G RS 1915+105 in a soft phase was analyzed in 
detail by lUeda et all ((2009); some simple photoioniza- 
tion modeling techniques were applied to describe the 
disk wind that was detected. The broadband spectral 
fits detailed in that work were used to calculate a bolo- 
metric luminosity for this observation (L = 7.5 ± 0.8 x 
10 38 ergs s" 1 ). As wit h the HETGS obse rvation of GRO 
J1655-40 reported in IMiller et al.1 (|2008[ ). this spectrum 
of GRS 1915+105 is particularly rich, and no single ve- 
locity can describe all of the ions observed. In this work, 
we adopt a value of v = 300 km s _1 because it is consis- 
tent with many ions and achieves a balance between the 
most and least ionized component s of the flow . The line 
properties and analysis reported in lUeda et al.l ([2009D are 
consistent with ft = 2.5 and log(£) = 4.3 (errors are not 
reported). We have used these values in calculating the 
k inetic power of t he win d in this observation. 

INeilsen fe Led (|2009) treat a number of Chan- 
riro/H ETGS observations of GRS 1915+105. lUeda et al.l 
(2009) focused on t he observation called "SI" in the 
INeilsen fe Lee (2009j) scheme; it is the lowest luminos- 
ity of five "soft" observations considered in their work. 
In each soft observation, an ionized X-ray disk wind is 
detected at high significance. To understand the range of 
wind properties in this important source, then, we have 
included observations S2-S4 in this analysis. 

Reduced and calibrated spectral and response files for 
each observation were obtained through the Chandra 
"tgcat" facility (for details of the observations, please 
consult Neil sen fe Led [20091) . The combined first-order 
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HEG spectrum from each observation was fit in the 2.3- 
9.0 keV band with a simple phenomenological model 
consisting of disk blackbody and power-law components. 
The lower energy bound was set by the high column den- 
sity along the line of sight to GRS 1915+105 (Njj was 
fixed at 5.0 x 10 22 cm 2 in each case); the upper bound 
was set by the likely calibration residuals in the high en- 
ergy portion of the spectra. This continuum model is not 
unique, but it allows for a simple and accurate charac- 
terization of the flux in the observed energy band, and 
gives reasonable basis for extrapolating to the 0.5-10.0 
keV band. 

We then calculated and applied grids of XSTAR mod- 
els. A customized grid was made for each observation, 
using the observed spectral continuum and the unab- 
sorbed 0.5-10.0 keV luminosity as the spectral input 
(the power-law was truncated below 1 keV to prevent 
unphysical results). Aga in, additiona l details of this pro- 
cedure can be found in IMiller et all (|2006al I2008D and 
IKing et al.l (l2012al lbh . In all turbulent velocity 

of 500 km /s, an iron abund ance of 2.0 times the solar 
value (e.g. lUeda et all [2009;. a density of log (n) = 12.0, 
and a covering factor of Q = 2.5 were assumed. In esti- 
mating the kinetic power in the wind for o bservations S2- 
S5, we used the velocity shifts reported bv lNeilsen fc Led 
(2009) , and the ionization parameters measured through 
direct fits to the Chandra spectra with the XSTAR mod- 
els. The values for all relevant parameters are given in 
Table H 

4U 1630—47: Last, we considered an archival Chan- 
dra HETG observation of the black hole candidate 4U 
1630-47. Observation 4568 started on 2005 August 9 at 
20:16:02 (UT), with a duration of 50 ksec. The instru- 
mental configurati on and modes used were the same as 
those described in [Miller et~al1 (|2008[ ). We again down- 
loaded the calibrated first-order gratings spectra and re- 
sponses using the Chandra "tgcat" facility, and generated 
a combined first-order HEG spectrum. 

Fits to the continuum with an absorbed disk blackbody 
plus power-law model give a high column density (N# = 
7.8xl0 21 cm" 2 ), a fairly hot disk (kT = 1.36±0.01 keV), 
and a power-law index of T = 2.00 ± 0.02. (The power- 
law index was checked by making fits to a simultaneous 
RXTEPCA spectrum over the 3-30 keV band.) Assum- 
ing a distance of 8.5 kpc, this continuum model gives a 
bolometric luminosity of 2.2 ± 0.2 x 10 38 ergs s — 1 . We 
further assume a mass of 10 M Q for 4U 1630—47 in this 
work. 

Again as per the procedure in IMiller et al.l (|2008l ). this 
spectrum was used to illuminate gas in a grid of XSTAR 
models. Solar abundances were assumed, and a turbu- 
lent velocity of 500 km s" 1 was found to give the best 
fits to the data. As with H 1743—322, the parameters 
of the binary system are not well known, and a covering 
factor of ft = 2tt was assumed in the generating the pho- 
toionization models. A fiducial density of log(n) = 12 
was also assumed in generating the models. As per the 
high luminosity observation of GRS 1915+105, only an 
H-like Fe XXVI line was detected, immediately indicat- 
ing a high ionization. Direct fits with the XSTAR grid 
give an ionization of log(£) = 4.9 ± 0.4, and a blue-shift 
of v = 300 ± 200 km s _1 . 



2.4. Jet Power 

In collecting a jet sample, we also wanted to create 
a set of uniform standards and conditions that would 
ensure high quality and rigorous results just as we had 
done for the wind sample. Jets are found in the ra- 
dio as a result of synchrotron radiation, so it is tempt- 
ing to utilize the radiative portion of the energy as 
an e stimate for the jet power (e.g., iMerloni fc Heinzl 
2007). However, the majority of the energy carried 
off by the je t s app e ars to be mechanical, n ot radiative 
dHeinz et al.l 120021: iDi Matteo et afl 120031: iGallo et al.l 
2005HAllen et alj|2006t i Tavlor et al.ll2006irWe note t hat 



Merloni fc Heinzl ([20071 ) and iCavagnolo et al.l ()2010D do 



find a relation between radio emission and mechanical 
power in jets, but a direct determination of the power is 
preferr ed over a proxy s uch a s radio luminosity. In ad- 
dition, ICavagnolo et al.l ([2010D report the jet power rela- 
tion to have a scatter of 0.70 dex and to be calibrated to 
high luminosity sources. It is unclear if extrapolation to 
lower luminosity sources is applicable. Likewise, the ra- 
dio l uminosities of jets are al so subject to Doppler boost- 
ing (jUrrv fc Padovanill991| ). which can be difficult to 
de- project, since the intrinsic spectrum must be known. 
We therefore restrict ourselves to only the most tentative 
of comparisons to radio luminosity via the fundamental 
plane of accretion onto black holes (IMerloni et afl 120031 : 
iFalcke et aLll2004t iGiiltekin et al.ll2009aft . This plane re- 
lates the radio luminosity of SMBH to the accretion rate 
(via X-ray luminosity) and mass of the black hole. See 
Section [4.51 for further discussion. 

For a more stringent comparison with our wind sam- 
ple, we require that the jet power be a direct estimate of 
the mechanical energy, not an indirect estimate using the 
radiative luminosity as a proxy for jet power. One way 
of quantifying the amount of mechanical power released 
via jets is to look at the volume they carve out in the 
form of "cavities" or "bubbles" . These cavities are seen 
in both the radio and X-ray wavelengths. The energy 
(Ej et ) is then estimated to be the sum of the internal 
energy and the PdV work done to inflate the bubble. 
The time needed to carve out such a region is estimated 
at i a9e =R/c s , where t age is the age of the bubble, R 
is the distance from the black hole to the center of the 
cavity and c s is the sound speed, typically estimated us- 
ing X-ray observations. There are a number of different 
ways to estimate the age of the bubble, but using sound 
speed and bubble radius is a fine approximation as long 
as t he bubble is still "attac hed" and is not buoyantly ris- 
ing ([Dunn fc F abian 2004]). Therefore, the power of the 
jet is Pj ot = Ej ct /t age . We note here that this estimate 
is a long term average and is not sensitive to discrete 
episodes of jet emission, which would serve to increase 
the jet power estimate. 

For this work, we d r aw di rectly from the sample de- 
scribed in lAllen et al.l ([20061 ) , who use this prescription 
to analyze a sample of nine elliptical galaxies that display 
such X-ray and radio cavities. They also estimated the 
Bondi mass accretion rates by constructing radial tem- 
perature profiles close to the black hole from Chandra X- 
ray observations. As the lAllen et al.l ([2006[ ) is a study of 
elliptical galaxies at low accretion rates, spherical accre- 
tion, i.e. Bondi accretion, is assumed with an efficiency 
conversion between mass accretion rate and luminosity 
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(77 = 0.1 as given in Allen et al. 2006). Conversely, in the 
wind sample, the accretion rates are typically higher at a 
few percent of Eddington, and the accretion is assumed 
to be through a standard thin disk. We u se both the 
j et po wer and Bondi luminosities reported in I Allen et al.l 

(2006) . We used the Bondi luminosities instead of the 
bolometric luminosities because each is an appropriate 
estimate for the mass accretion r ate in these particular 
systems. iMerloni fe Heinj (120071) rep ort analysis of the 
same nine AGN as lAllen et all (|20M ) as well' as six ad- 
ditional sources. We do not include these extra sources 
in our sample as the temperature inside the Bondi ra- 
dius was not measured directly but extrapolated from 
much further, outside regions. Howeve r, as will be shown 
in Section 14.21 IMerloni fc Heind (120071) repo rt consistent 
analysis and results with lAllen et al.l (|2006| ) . 

At the low mass end, it is much more difficult to use 
the same methods to estimate the jet power from bub- 
bles and cavities. This is because most of the black hole 
candidates are not in regions with dense gas, making the 
bubbles hard to observe. The few that are embedded in 
dense clouds happen to also be in star forming regions. 
This means that the observed cavities can be carved out 
not only by their jets, but also by their high mass com- 
panion star's winds or even the supernovae associated 
with the black holes themselves. One such candidate is 
BHB, Cygnus X-l. This is a stellar mass black hole with 
an associa ted radio bubble that is thought to b e created 
by its jet (|Gallo et al.ll2005t iRussell et aLll2007h . This is 
evidenced by the fact that the long axis of the bubble is 
aligned with the jet axis. However, there is no counter- 
jet seen, and the b ubble may also be a ssociated with a 
s upernova remnant (IRussell et al.ll2007fl . 

IGallo et al.1 (|2005| ) and lRussell et al.1 (|2007| ) have both 
used the observed cavity to estimat e the power of th e 
jet, employing similar techniques as A llen et al.l {2006). 
We note that in making dens it y esti mate s of the emit- 
ting r egion, both IGallo et "ail (|2005j) and Russ ell et al.l 

(2007) assu me the emission is B remsstrahlung radiation. 
However, as lMarti et al.l (|1996| ) show, the radio emission 
in the limb brightened areas have steep spectra, which 
impli es the regions are n on-thermal in nature. More- 
over, IRussell et al.l (|2007[ ) find that the emitting loop 
is not visible in the R band, indicating that the emis- 
sion detected with an Ha filter, is in fact Ha rather 
than Bremsstrahlung emission, as assumed by Gallo et 
al. (2005). 

One can use the observed Ha flux to estimate the av- 
erage density o f ionized ga s (abou t 6 cm -3 ) and follow 
the method of IGallo et al.l (|2005| ) to determine the jet 
power, Pjot = 10 34 — 10 38 ergs s _1 . As is obvious by the 
four orders of magnitude, there are large uncertainties 
that go into this calculation. Clumpiness of the emitting 
gas would give an overestimate of the average density, 
and if the ionized gas is indeed produced by a 100 km 
s _1 shock, it occupies a very thin sheet compared to the 
apparent size of the emitting region. Second, the neu- 
tral fraction in th e emitting reg ion, taken to be 98% by 
IGallo et all (|2005D and zero by IRussell et all (|2007l ). is 
not well known. Finally, when shocked gas has cooled 
to the point that Ha emission is efficient, its pressure is 
probably dominated by the magnetic field, so the sound 
speed should be replaced by the fast mode speed, which 



is several times larger. 

Instead, if we comb ine the intensity measurements 
of IRussell et al.l (120071) with the shock wave models of 
iRavmond et al.l (|1988| ) and the theory of int erstellar bub- 
bles blown by a continu ous energy input (|Castor et al.1 
ll975l : IWeaver et al.lll977l ) , we can get a t ighter constraint 
on the power estimate. IRussell et al.l (|2007t ) measured 
an intensity in the [O III] AA5007, 4959 lines of approxi- 
mately 1.5 x 10~ 15 ergs cm~ 2 s _1 per square arcsec after 
correction for extinction, in a 2' section of a cut through 
the NE part of the shell. They also measured [O III] to 
Ha + [N II] ratios that indicate shock speeds of 90-200 
km s _1 (their figure 7). Shock waves in that range pro- 
duce 0.87±0.2 photons in the [O III] lines per H atom 
that passes through the shock (|Ravmond et al.| [l988). 
The shell is limb-brightened, and comparison of the 2' 
thickness with the 5' radius indicates an enhancement 
factor of 2.4. Thus 

1.7 x 10 7 = 2.4 * 0.87 * n iy (47r)photons cm" 2 s" 1 sr" 1 

w 

where no is the pre-shock density and V s is the shock 
speed in units of cm s -1 . Thus uqV s = 1.0 x 10 8 . The 
expression for the radius of a wind-blown bubble in the 
intermediate stage (when the shock is radiative as in the 
Cygnus X-l bubble; equation 21 of iWeaver et al.l 119771) 
can be converted to 



L 37 = 7.7 x 10- 8 n V;V 1 2 00 i? 2 (2) 

where L37 is the jet luminosity in units of 10 37 erg s _1 , 
UqVs is in units of cm s~ l , T^oo is the shock speed in units 
of 100 km s _1 and R5 is the bubble radius in units of 5 
pc. For a shock speed of 90-200 km s _1 and an average 
bubble radius of 4 pc, this impies a jet luminosity of 
4 - 20 x 10 37 erg s _1 . 

There are two important caveats to keep in mind, both 
of which could lead to a severe overestimate of the jet lu- 
minosity. First, O stars in the region, including the com- 
panion of Cyg X-l itself, might contribute ionizing flux 
that we are assuming to come from the bubble shock. 
Second, the shell could be a result of the explosion that 
created the Cyg X-l black hole, rather than the jet. The 
alignment of the shell with the jet direction suggests, 
however, that the jet plays a significant role in the en- 
ergetic of the bubble. In addition, there is some uncer- 
tainty involved with the reddening correction. In view of 
the uncertainties in the relevant quantities, in the best 
method of estimating the power required to inflate the 
bubble, and the origin of the bubble itself, we use the full 
range of power estimates noted above. We note there are 
tighter constraints if the bubble is only associated with 
inflation by the jet and reddening is unimportant. 

2.5. Ultra-Fast Outflows 

Finally, we also defined a smaller subsample of winds 
that are moving faster than 3,000 km s _1 (0.01c) rela- 
tive to the systemic velocity. For this sample, we relax 
our standards for the AGN and include both Suzaku and 
Chandra imaging spectrometers. As we only discuss four 
examples of these outflows, (one BHB, one nearby quasar 
and two gravitationally lensed, higher redshift quasars), 
this additional data set is only meant to be illustrative 
not exhaustive. This sample still requires that the ab- 
sorption features be at least 3cr significance. 
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The first of these ultra fast o utflows is the BHB 
J17091-3624. IKing et al.l (|2012b| ) discus two Chandra 
HETG observations, one of which has clear absorption 
features above 6.9 keV. Using XSTAR photoionization 
grids, they model these features self-consistently and find 
an ionization of the absorbing gas to be log£ = 3.3^q'^, 

moving at v ont = 9600^00 km s _1 . We also use the 
second component at a slightly higher ionization state 
log£ = 3.9±g;| and velocity, v out = 15,400 ± 400 km 
s -1 . These are the fastest outflows observed from a BHB 
candidate, and they bear re semblance to some of the 
outflows seen in qua sars (e.g., IChartas et aUl2"002l [20071: 
ITombesi et aHl20ll. 

ITombesi et alJ ([20111 ) used the Suzaku X-ray Imag- 
ing Spectrometer (XIS) to observe 3C 111, an AGN at 
z = 0.0485. Of the three observations in their study, 
one showed evidence of an absorption line in the Fe K 
band. This observation had an exposure of 59 ksec. Us- 
ing XSTAR grids with turbulent velocity of 3000 km 
s _1 , this one feature was fit with an ionization parame- 
ter of log£ = 4.32 ± 0.12 and an outflowing velocity of 
v t = Q.1Q6 ±q .0 06c. 

Chartas et alJ ([2002) observed the quasar APM 



08279+5255 with the Chandra Advanced CCD Imaging 
Spectrometer (ACIS) and noted outflows in the absorp- 
tion spectra. By using lensed quasars, one is able to 
probe outflows which would otherwise be too faint to ob- 
serve. This quasar is at a redshift of z — 3.91. The 
spectra of APM 08279+5255 shows two features at 8.05 
keV and 9.79 keV in the rest frame of the host galaxy 
(|Chartas et all 12002ft . If these correspond to Fe XXV 
then the outflowing velocities would be 0.2c and 0.4c, re- 
spectively. We utilized b o th com ponents in this analysis. 

Finally, IChartas et al.l (|2007l ) described a gravitation- 
ally lensed quasar PG 1115+080, which is at a redshift 
of z = 1.72. The authors also used ACIS and notice ab- 
sorption features in the host Fe K band. PG 1115+080 
has prominent absorption features at rest frame ener- 
gies of 7.27 keV and 9.79 keV; both of which are used 
in our analysis. Associating these features with Fe XXV 
gives velocities of 0.09c and 0.40c. These features were 
not modeled with a photoionization model, but an as- 
sumed ionization of log£ = 3.5 is taken as an esti- 
mate for the ionizat ion parameter from XSTAR models 
(IChartas et al.ll2007ft. Bol ometric luminosities are taken 
from lChartas et al.l (j2007f l. while estimates of the ioniz- 
ing luminosities were taken as the rest frame X -ray lu- 
minosity from 0.2-10 keV from [DaLeFaD d20Q3). 

3. METHODS 

After acquiring the sample, we calculate the mass out- 
flow rate in the wind systems. This is done using simpli- 
fied, order of magnitude estimates based on the expres- 
sion for spherical wind. It is modified by both covering 
and filling factors to account for the fact that winds are 
not spherical outflows: 



ilpr vC v 



(3) 



where, is the covering factor ( < Vt < in), p is the 
mass density (p = pm p n e ), p is the mean atomic weight 
assumed to be p — 1.23, m p is the mass of a proton, 
n e is the electron density, r is the radius from the ion- 
izing source, v is the out-flowing velocity, and C v is the 



line of sight global filling factor. As the winds may be 
clumpy and filamentary, they are likely to have a small 
filling factor. This expectation is based on the observed 
varia bility of the abso rption lines (e.g.. | Crenshaw et alJ 
l2003tlElvis et al.ll2004tlRisariti et al.l2009ft aswellasden- 
sity diagnostics (e.g.. IKing et al.ir2012aj ). We note that 
variability can be due to both motion along our line-of- 
sight as well as the duty cycle of the wind. The short 
timescales of variability suggest that the variability is 
likely due to small filling factor and clouds moving across 
our line of sight rather then dissipation of the wind it- 
self. Further, there is an inconsistency in the literature 
as to the actual filling factors, as there are few direct 
constraints on this quantity. This factor should vary be- 
tween different sources as well as with ionization, £, but 
a range from 10~ 5 to 1 is seen in the literature across 
the full mass scale (e.g. iMiller et al l 120081: jMocz et alJ 
120111: iZhang et alJ 1201111 IKing et al H 2012allbT). Finally, 
as shown in iGiustini fe Progal ([2012D . the exact nature 
of the wind may depend on the ionization, velocity and 
density of the wind, and may be quite uncertain. There- 
fore, we do not assume a value for the filling factor for 
any of the measurements but leave the kinetic energy 
luminosity in terms of the filling factor. 

Equation [3] can be rewritten in terms of observable 
quantities from the X-ray absorption features using £ = 
Li on /n e r 2 , where £ is the ionization parameter and Li on 
is the ionizing luminosity between 1 - 1000 Ryd (1 Ryd 
= 13.6 eV). 

• _ pm p VlLi on vC v . 

jyio-at — 7 (4 ) 

For consistency, we use H, = 2ir for al l the AGN sources, 
based on findings by iRevnoldsl ([1997ft who found half of 
all Seyferts show evidence for warm absorbers. We have 
not assumed a filling factor, C v , but have reported our 
results of kinetic luminosity per filling factor. To con- 
vert the mass outflow rate to the kinetic energy carried 
away by the warm absorbing winds, i.e., power or kinetic 
luminosity, we use the following relation, 



1 • 2 

^wind = -M out V — 



pnfipQLionV^ C v 
2£ • 



(5) 



The total kinetic luminosity is the amount of mechani- 
cal energy that is carried away by the wind. It is im- 
portant to understand this in the context of the total 
escaping energy of the black hole and accretion disk. 
This can be done via comparison to the radiative power 
released, i.e., the bolometric luminosity. Further, the 
bolometric luminosity is often considered a proxy for the 
mass accretion rate by assuming an efficiency conversion, 
£bo1 = p Mc 2 . The efficiency , 77, is usually assumed to be 
10% (e.g.. lAllen et al.|[2006t IVasudevan fe Fabianl [2001 
iFabian et al.ll2009ft which is consistent with the Soltan's 
argument "(ISoltanlll982T) . However, in reality it is likely 
to vary between sources and with Eddington fraction. 
In our study, the bolometric luminosities for the AGN 
are taken from broad-band spectral energy distr i butio n 
(SED) fitting performed by IVasudevan fe Fabianl (|2009ft . 
In the few instances that the AGN lack a bolometric 
luminosity, we used the co nversion Lb i ~ 20L2_io/ceV 
(jVasudevan fe Fabian 2009). For the stellar-mass black 
holes, their SED peaks in the X-ray. Therefore, the bolo- 
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metric luminosities are taken from the X-ray observations 
as the luminosity between 0.5-10.0 keV. Where values 
in the literature were quoted for different energy bands, 
the luminosity was converted to the 0.5-10.0 keV band 
for consistency by extrapolating the given models within 
Xspec. We also note that any uncertainty in distance, 
which could effect the bolometric and ionizing luminosi- 
ties as well as estimates of the outflowing velocity, are 
small as compared to the uncertainties in these measured 
quantities. 

The jet power is calculated using the energy estimates 
of radio and X-ray cavities and age of the bubble as de- 
scribed in Section [HH Pj e t = Ej ei /t age . We again note 
that this estimate is a long term average and short, dis- 
crete episodes of jet emission would increase the jet power 
estimates. 

4. ANALYSIS AND RESULTS 

4.1. Bolometric Luminosity versus Wind Power 

After acquiring a sample of BHB and AGN with esti- 
mates for the wind power, we begin to analyze how LboI 
relates to the total kinetic outflowing power in each sys- 
tem. In this initial analysis, we only include the lower 
velocity winds (not the jets, which we consider in Sec- 
tion [4T2l or the ultra-fast outflows, which we consider in 
Section 14. 4|) . Figure Q] shows the distribution of veloci- 
ties and ionization parameters that are included in the 
kinetic power of the winds. Figure [2] shows the kinetic 
wind luminosity as compared to the source bolometric 
luminosity for individual components in each observa- 
tion, while Figure [3] shows the same plot but for the 
total kinetic wind luminosity for each observation. The 
stellar-mass black holes cluster at the lower luminosities, 
while the SMBH are found at the higher luminosities as 
expected. A positive correlation is apparent in the data 
set. 

In Figure [2] (as well as the following figures), we plot 
the kinetic luminosity per filling factor, C v , versus the 
bolometric luminosity. Further, Figure |3] shows the total 
kinetic power for each observation, which uses the sum 
of the individual components plotted in Figure [2] We 
begin our analysis with the total kinetic power for each 
observation. Initially, we assume that there is a common 
relation between both the AGN and BHB. In Section l4~5l 
we relax this assumption and characterize the two groups 
separately. 

To first characterize the trend given in Figure [3J we 
utilize two correlation tests: a Spearman's rank test and 
a Kendall's tk test. We find a Spearman's rank coeffi- 
cient ps = 0.89, with a null hypothesis probability (i.e., 
no correlation) at p = 1.8 x 10 -8 , indicating a strong and 
positive correlation. The value for Kendall's coefficient is 
tk — 0.72 with p — 1.4 x 10~ 6 , also indicating a strong, 
positive correlation. 

Following this, we assume that the data can be de- 
scribed by the linear relation, 

y = ax + (5 (6) 

where, y = \og(L wiadA2 /C v ), x = log(i Bo l,42), and the 
subscript "42" denotes the units 10 42 erg s" 1 . We then 
minimize the function, 

2 = (?/»-/3-az») 2 (7] 



to estimate a an d (3 (e.g., iPress et all 119921 : 
iTremaine et al.l 120021 ) . Here o~ y i and a X i are the 
errors associated with the kinetic wind luminosity and 
the bolometric luminosity, respectively. The quantity o~o 
is the intrinsic scatter in the relation and is determined 
by ensuring the reduced \ 2 is close to unity. We obtain 
a = 1.58 ± 0.07, /3 = -3.19 ± 0.19 and a a = 0.68, such 
that, 

log (£ w ind,42/C'u) = 

(1.58 ± 0.07) log (L Bo L4 2 ) ~ (3.19 ± 0.19). (8) 

These parameters are listed in Table [3J The large <To 
implies that the intrinsic scatter in these measurements 
is dominant over the measurement errors. We can expect 
a high intrinsic scatter due to the high variability of each 
of these sources, especially because the observations used 
to derive the AGN bolometric and wind luminosities of 
individual sources were not made simultaneously. Envi- 
ronment may also play a large role in driving this scatter, 
evidenced by the fact that the larger scatter is associ- 
ated with the SMBH measurements, which are located 
in dense groups and clusters to open field environments. 
This scatter may also be attributed to the bolometric 
correction applied to the X-ray luminosities of the AGN 
and not to the BHB which peak in luminosity in the 
X-ray band. 

Finally, the scatter may be due to the inclusion of a 
range of ionization parameters, especially in the AGN. 
There appears to be a stratification in the power which 
may depend on the ionization parameters. This is shown 
in Figure [2] as the low ionization components tend to 
have higher powers as compared to the higher ioniza- 
tion parameters. Therefore, we also examine only the 
high ionization, i.e. log£ > 2, components with an aim 
to better compare the same sample in both AGN and 
BHB. We fit the individual components with high ioniza- 
tion and find the slope flattens to a log?>2 = 1.42 ± 0.06, 
/3 iogS>2 = _(3.734-0.14) and a scatter ofa° ei>2 = 0.57. 
These results are given in Table [3] and shown in Figure 

m 

4.2. Jet Power 

We next include the relation between the Bondi lumi- 
nosity and jet power as estimated via the radio bubbles 
seen in elliptical galaxies and Cygnus X-l (See Figure [4] 
and Table [2]) . Again, the Bondi luminosity is used for 
the ellipticals instead of bolometric luminosity as in the 
wind sample, but can be thought of as a proxy for mass 
accretion rate just as the bolometric luminosity is at high 
accretion rates. For the jet sample, a high degree of cor- 
relation between the Bondi lumi nosity and the jet power 
is indicated, as initially noted bv I Allen et al.l (|2006l ). We 
find a Spearman's rank coefficient of ps = 0.95 with 
probability p = 2.3 x 10 -5 , and a Kendall's tk coefficient 
of t k = 0.87 with probability p = 4.9 x 10~ 4 . Clearly, 
there is a positive correlation for this data set as well. 

As a correlation is quite apparent in this jet sam- 
ple, we fit the data using the same technique as was 
used for Sections O We find a a jet = 1.18 ± 0.24, 
/3 jet = —(0.96 ± 0.43), and a intrinsic scatter consistent 
with zero ((Tq 6 * = 0). These parameters are listed in Ta- 
ble [3] for comparison with the wind parameters. This 
relation is also shown in red in Figure [4] The reduced 
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X 2 is quite small at \ 2 jv = 0.11. This is a result of 
the large uncertainty estimates on the Bondi luminosi- 
ties. The error in Bon di luminosity is e stimated using 
the uncertainty given in I Allen et al. (2006) as well as the 
uncertai nty derived from the scat ter in the M — a relation 
given by iGiiltekin et al.l ([20"0 9bl . These two uncertain- 
ties are added in quadrature, resulting in an uncertainty 
of a pproximately 0.62 dex. 

As lAllen et atl ((20061 ) and lMerloni fc Heinzl ([20071 ) per- 
form similar fits to exclusively the elliptical galaxies, we 
next exclude Cygnus X-l from the fit and find a? et = 
1.34 ± 0.50, /3 jet = -(0.80 ± 0.82) and an intrinsic scat- 
ter also consistent with zero. See the orange line in 
Fi gure |H Our analys is is able to reproduce the results 
by Allen et al.l (120061) . who first published this sample. 
lAllen et al.1 (j2006[ ) find B = 0.77 ± 0.20, which is equiv- 
alent to a jet — 1.30 ± 0.34 in our nomenclature . Our 
results are also consistent with iMerloni fc Heinzl ([20071 ) 
who find a slope of a? et = 1.6^*3 when correlating 

Lifm/iEdd to Lsondil '-^Edd- 

We note here that Cygnus X-l nominally lies one or 
two orders of magnitude above the elliptical jet rela- 
tion but is consistent within its large uncertainty. This 
begs the question whether the radio bubble seen around 
Cygnus X-l is truly related to the black hole jet or in 
fact a chance alignment. Cygnus X-l is located in a 
fairly active star forming region where mas sive young 
stars may be responsible for such a structure (|Reid et al.l 
120111 ). In fact, there have been X-ray winds associated 
with the Cygnus X-l system, whether from the accre- 
tion disk or the companion O star is unclear. These may 
also have an effect on inflating the bubble, which would 
bring the power estimate down. In addition, as previ- 
ously discussed, we do not see a counter bubble from the 
presumed counter-jet. 

If we include Cygnus X-l in the jet relation, then the 
slope of the jet relation is inconsistent with the wind re- 
lation by only 1.6<r. However, when we exclude Cygnus 
X-l ( as it is plausibly as sociated with a supernova rem- 
nant iRussell et al.ll2007l) the slopes of the jet and wind 
relations are consistent within errors. Although the nor- 
malizations are different, a common slope might imply 
a shared driving mechanism. In Figure |4j it is apparent 
that the jet and wind power normalizations are within 
a few orders of magnitude, especially at high luminosity. 
However, correcting for the filamentary and geometric 
structure of the winds via the filling factor will decrease 
the wind power normalization by 3 to 4 orders of mag- 
nitudes, demonstrating the greater efficiency of the jet 
power. 

4.3. Spectral State Dependence 

It is also interesting to compare how jet and wind 
power scale in terms of Eddington fraction to examine 
the accretion rate dependence. Figure [5] shows both the 
wind and jet power per Eddington luminosity as com- 
pared to their Eddington fraction (i.e., bolometric lumi- 
nosity or Bondi luminosity per Eddington luminosity). 
The AGN jet power is denoted in red, Cygnus X-l jet 
power is denoted in orange, the AGN winds are denoted 
in black, and the BHB winds are denoted in b lue. The 
solid lines are taken from lChurazov et al.l (|2005l ) and de- 
scribe the outflow mechanical power (thick line) and ra- 



diative power (thin line) for AGN. They postulate that 
AGN should follow a similar evolution to their stellar- 
mass counterparts, in that they should have a strong jet 
dominated phase at low accretion rates, and little-to-no 
jet production at high accretion rates. This allowed them 
to present a model for AGN feedback and co-evolution 
with the host galaxies as a function of mass accretion 
rate. We assume an efficiency of ?y = 0.1 to compare our 
bolometric luminosity to their mass accretion rate. 

A division in the outflow power is seen at approxi- 
mately 10 _2 LEdd in our data set. See Figure [5] Below 
this Eddington fraction, jets dominate and increase in 
power with bolometric luminosity. At higher Eddington 
fractions, the wind power dominates but decreases with 
Eddington fraction. We note that there is an observed 
dichotomy between the type of outflows seen in stellar 
mass black holes and their X-ray spectral state already 
described in the literature. In particular, winds are found 
in the "high/soft" state, i.e., high mass accretion rate and 
Eddington fractions (> 10~ 2 LEdd), and radio jets are ob- 
served in the "low/hard" stat e, i.e. low mass accretion 
rates and Edd i ngton fractions (|Miller et al.ll2006al 120081: 
iNeilsen fc Led 12009ft . Although our sample is not ex- 
haustive, we do see illustration of a similar trend in the 
AGN sample where jets persist at low Eddington frac- 
tions and winds persist at high Eddington fractions. We 
note that in one AGN source, NGC 4051, there is ev i- 
dence for simultaneous winds and jets ([King et al.ll2011l) . 
However, the winds may very well dominate at this high 
Eddington fraction in NGC 4051. Regardless of the par- 
ticular outflow seen, Figure [5] allows for the prediction of 
the outflow power as a function of Eddington fraction. 
This will be vital to simulations of AGN feedback and co- 
evolution; both for matching predictions to observations 
as well as implementing sub-grid physics in cosmological 
simulations. 

4.4. Ultra-Fast Outflows 

There are a few wind sources that lie well above the 
wind relation in Figure [3] and much closer to the jet re- 
lation in Figure 01 This is primarily because of their 
high velocities, i.e., v > 0.01c, as L W md °c v 3 /£, which 
increases their power to lie near the jet relation. In Fig- 
ure [SI we include four ultra-fast outflows mentioned in 
Section [2751 These are denoted in the black squares. The 
upper squares assume a global filling factor of unity. For 
comparison, we include estimates of the wind luminos- 
ity if the filling factor were as low as C v = 10~ 4 con- 
nected by a dashed line to the original higher estimate. 
Such a small filling factor is reasonable as it is consis- 
tent w ith their potential transient na ture (Tombcs i et al.l 
2011) as well as density diagnostics ([King et al.M 2012a) . 
We note that variability of these sources can be at- 
tributed to both movement across our line-of-sight as 
well as duty cycle. However, the timescales for vari- 
ability are short compared to dissipative timescales, and 
are gener ally ascribed to fi l ling factor and not duty cy- 
cle (e.g.. lElyis et all [2001 ICrenshaw fc Kraemerl 120121: 
IRisaliti et al.1120091 ). Regardless of the filling factor, these 
high velocity outflows tend to be much more efficient 
at their given bolometric luminosity as compared to the 
other winds. Therefore the ultra- fast outflows may re- 
semble the jet relation, which is plotted in Figure [6] As 
the ultra-fast outflow power approaches that of jets, it 
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suggests that we are seeing the transition from winds as 
they are being accelerated into jets. 

4.5. Distance and Mass Dependence Diagnostics 

When examining broad relationships, it is important 
to be wary of a rising trend with a slope of unity; this 
can indicate that a relation that is driven by a mutual 
dependence of a third, shared parameter. For example, 
in the wind relation, both the bolometric and ionizing 
luminosity are both proportional to the square of the 
distance, which may have an influence in driving the ob- 
served trend. However, we find evidence to the contrary. 
Not only is the observed slope in the wind relation is 
greater th an unity, but a partial c orrelation test, de- 
scribed by lAkritas fc Siebertl (|1996 ), gives a low prob- 
ability of p ~ 0.035 that the wind relation is driven by a 
mutual dependence on distance. This estimate is derived 
from Kendall's partial tk, p , which gives the Kendall's r 
holding the third parameter, distance, constant. The val- 
ues we find for the wind relation are tk, p — 0.270, with 
an estimated variance, <jk — 0.128. Further, the rising 
trend in the winds is therefore dominated by the veloc- 
ity and ionization as L W md oc u 3 /£, both of which are 
directly and independently constrained by observations. 
We also note that the jet sample is not driven by dis- 
tance. This sample has an even smaller probability of 
p < 10~ 6 that distance is needed as a third parameter. 
This is derived from a Kendall's partial correlation test 
where r^ p^O.86 and o\r-=0.13. 

Similarly, we tested whether the relations are driven 
by the mass of the black hole. A partial correlation test 
of the wind sample using mass as a third variable gives 
a probability of p = 0.034 that the relation is driven 
by a mutual dependence on mass. This probability was 
derived from a Kendall's partial correlation test where 
we found tk, p is 0.436 with ok — 0.206. The jet sam- 
ple is even less dependent on mass with a probability of 
p = 1.8 X 10~ 3 (t k , p = 0.73 and a K = 0.24), likely the 
result of the sample including only 1 BHB and massive 
ellipticals (M ~ 10 9 M BH )- 

Although the probability in the wind sample is small, 
it does not rule mass out as a third parameter at a 3cr 
level like the jet sample does. The nature of these winds 
across such a large mass scale has not been studied be- 
fore. Consequently, we explore the potential for mass 
dependence in this data set. 

We separately fit a linear relation to both the BHB 
and AGN data sets. We used the same linear regression 
as in Sections 14.11 & 14.21 minimizing the x 2 for the best 
fit parameters. We find the BHB sample had the best fit 
parameters of a BHB = 0.91±0.31, /3 BHB = -(5.58±1.68) 
and an intrinsic scatter consistent with zero (cr BHB =0). 
The AGN sample had the best fit parameters of a AGN = 
0.63 ±0.30 and /3 AGN = -(1.24 ±0.63) with an intrinsic 
scatter of cr AGN = 0.58. Figure depicts the best fit 
relation for the BHB and AGN, and Table [3] lists these 
parameters for comparison with previous results. 

The slopes of the two individual fits are consistent with 
each other, although inconsistent with the initial fit to 
the wind sample at the 2.4cr and 3.1cr level for the BHB 
and AGN samples, respectively. However, in order to 
evaluate whether these parameters are truly inconsistent 
with the initial fit, we used a bootstrap method to re- 



sample the data and estimate the number of trials we 
would expect with slope a > 1.58. Using N=10 4 trials, 
we found that in the BHB sample, 3.1% of the trials gave 
a slope a > 1.58. Similarly, in the AGN sample, 1.7% of 
the trials gave a slope a > 1.58. Although these parame- 
ters are formally inconsistent with the initial fit using the 
1 cr error bars, we can not rule them out at more than a 
98.3% confidence level. In addition, by fitting the data 
separately, we are introducing an additional three free 
parameters and thereby doubling the parameter space. 
Because we can not directly compare the Ax 2 , as the 
\ 2 jv is fit to be unity, we used a Bayesian analysis to 
determine which model better describes the data. 

The advantage to using Bayesian statistics is that it al- 
lows us to compare two different models of the same data 
without a reduced x 2 an d without the same number of 
degrees of freedom. This is done via a Bayesian odds 
ratio, which compares the likelihood of each model over 
the entire parameter space. In our analysis, we assume 
a uniform prior distribution in slope (a € [—10, 10]), in- 
tercept (/3 e [-12,2]), and scatter (a € [0,3]). We find 
that when comparing the likelihood of the single fit, £i, 
to the two individual linear fits, £2, the odds ratio was 
01.2 = g- = 0.17. This means that the two linear fits 
are slightly favored over the single linear fit, suggesting 
that mass may have a role in this relation. In Section \5. 31 
we suggest further observations that would also help to 
distinguish between the two models. 

We next fit a plane to the data, which included mass 
as an additional parameter. A similar linear regression 
to that used in Section B~T1 is used to fit a plane to the 
entire wind data set. The plane is described by, 

Z = a p X + (3 p Y + lp (9) 

where X = log(L Bo i), Y = log(M BH ), Z = 
log(L w ; n d/C 1 ,), and 7 P is the intercept. In order to find 
the best fit parameters, we minimized the function, 

2 = ( z i ~lv~ PyYj - a P Xi) 2 , , 

i=l a p"X,i + Pp a Y.i + "p,0 

where ax,i is the log(LBoi) scatter, ay,i is the log(AfeH) 
scatter, and <t Pj o is the intrinsic scatter of t he plane in the 
Z dire ctio n. This is further discu ssed in iMerloni et al.l 
(2003]) and lGultekin etafl (|2009aD 

We found the data set is best fit by the parameters are 
a p = 0.2 ± 0.4, P P = 1.2 ± 0.3 and j p = 24.5 ± 0.2 with 
an intrinsic scatter of a p .o = 0.68. Figure [5] shows the 
best fit plane. 

Ideally, we would like to compare this plane to a plane 
that describes the jet cavities while including mass. How- 
ever, the jet cavities are dominated by ellipitical galax- 
ies with masses of approximately Mbh ~ 10 9 M Q , and 
as demonstrated by the partial correlation test, the data 
show a relation that is independent of mass at the 99.82% 
confidence level. Although mass may be an important 
variable, the data may not span a wide enough parame- 
ter space to deduce its effects. However, previous studies 
that have used radio luminosity instead of X-ray cavi- 
ties to study jet characteristics, hav e shown a mass de- 
pendence (e.g.. IMerloni et al.l 120031: iFalcke et al] l2004t 
Gtilt ekin et all l2009al ). As we noted in Section 12.41 us- 
ing radio luminosity to study jet properties involves some 
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uncertainties whe n converting between ra dio flux density 
to jet power (e.g. Merloni & Heinz 2007), and should be 
treated with caution. 

Therefore, we proceed with only a tentative compari- 
son of our relation to the fundamental plane of accretion 
onto black holes. The wind plane parameters, a and f3, 
are formally consi stent with the fundame ntal plane pa- 
rameters given bv lGultekin et al.l (|2009aD . although the 
overall normalizations differ. See Figure [9l The plane 
given bv lGultekin et all (|2009al ) relates the accretion rate 
of low-luminosity black holes via the X-ray luminosity 
and the mass of the central black hole to the radio lumi- 
nosity of the compact radio source in the host galaxy. We 
note that we do not include a conversion between radio 
luminosity and jet power in this comparison, which may 
be necessary for com parison with our wind sample. In 
addition, the work bv lGultekin et all (|2009af ) uses X-ray 
luminosity and does not include the bolometric correc- 
tion factors that we have used in our work. Both of these 
caveats should be examined in the future to better un- 
derstand the connection and consistency between these 
two relations. Finally, we note that our wind coefficient 
describing the bolometric lu minosity, a n , i s cons istent 
with the coefficient given bv iMerloni et al.l (|2003| ). and 
our mass coefficient , f3 p , is only inconsistent with the 
Merloni et al. | (|2003[) coefficient at the 1.3 a level. Again, 
Merloni et al.l (|2003l) describe a similar plane between the 
mass, X-ray luminosity and compact radio emission in 
both stellar- and supermassive black holes. 

5. DISCUSSION 

We have compiled samples of both X-ray winds and rel- 
ativistic radio jets that span eight orders of magnitude in 
black hole mass. Each sample has uniform, rigorous se- 
lection criteria to ensure consistent comparisons between 
the various sources. In particular, we demand that the 
winds be detected through significantly blue-shifted ab- 
sorption features seen in the X-ray band, and the jets be 
seen as X-ray bubbles or cavities. By including only the 
X-ray winds and jets, we aim to probe the outflows as- 
sociated with the inner accretion disk. These flows may 
be driven by the accretion disk, and jetted outflows may 
also tap the spin of the black holes. 

For comparison, we also examine the bolometric or 
Bondi luminosity of each source. In doing so, we find a 
relation that describes the entire black hole X-ray winds 
sample as log(L w - md ,42/C v ) = (1.58 ±0.07) log(L B oi,42) - 
(3.19 ± 0.19) and jet sample as log(Lj ct 42) = (1.18 ± 
0.24)log(L B o„di,42) - (0.96 ± 0.43). If we exclude 
Cygnus X-l the relation becomes log(Lj ct ,42) = (1.34 ± 
0.50) log(i B ondi,42) - (0.80 ± 0.82). These relations sug- 
gest a common regulation scheme for winds and jets 
across the mass scale. We also find that when fit individ- 
ually, the BHB and AGN wind samples have shallower 
slopes of a BHB = 0.91 ± 0.31 and a AGN = 0.63 ± 0.30, 
which are also consistent with each other within errors. 
Although, the two wind fits are preferred slightly over the 
entire sample fit, a common slope between the BHB and 
AGN is required by the data regardless of the procedure 
used. 

5.1. Plausible Outflow Driving Mechanisms 

Examining winds specifically, thermal, radiative and 
magnetic mechanisms are viable methods of driving 



winds. However, it is not clear that these mechanisms 
would collectively drive this relation in the same way as 
required by the data. We now turn to whether radiative, 
thermal or magnetic processes can drive the observed X- 
ray wind correlation individually. 

First we examine whether radiation pressure, and more 
specifically UV line driving, has enough force to launch 
winds. If we assume this occurs when the force from 
the lines, Fu nes , exceeds that of gravity, F grav , i.e. 
Fimes > Fgravi A yields the following UV luminosity, 
Ljjy , criterion that LjryM (t) > L^dd (See equation 
in IProga fc Kallmanl|2"0 02). Here M(t) is the force 



multiplier (|Castor et al.lll975h . and L^dd is the Edding- 
ton luminosity. The force multiplier allows us to quan- 
tify the contribution of line driving in addition to elec- 
tron scattering and is a function of the optical depth 

p is the density, vth is the thermal 



+ ~ ™. I dv I 



velocity and 
(|Castor et al 



^ I is the velocity gradient along the flow 
1975). For most BHB, the strong X-ray 



radiation can highly ionize the gas, driving the M.(t) to 
1 at l ogg ks 2, and M{t) to 0.1 at log£ 3 (jProga et al.l 
|200d ). In addition, the BHB spectrum does not have 
a large relative contribution from the UV, due to the 
high disk temperatures . This also hinders line - driving 
of win ds consistent with lProga fc Kallmanl (|2002| ); |Progal 
( 2002). On the other hand, AGN spectra peak in the UV, 
and the AGN winds span a wider range of ionization pa- 
rameters, suggesting that at low ionization parameters, 
log£ < 2, line-driving may be importan t, which is also 
consistent with the work bv IProga] (|2002l) . This may par- 
tially account for the index of 1.58±0.07 in the initial fit 
that makes AGN more efficient wind producers. How- 
ever, it still remains to be seen what drives the higher 
ionization states found in a majority of the AGN listed 
here. 

Thermal pressure is another plausible driving mech- 
anism. Winds can be driven by thermal pressure if 
the temperature o f the gas is higher t h an the local 
escap e speed (e.g., IBegelman et al.l 1983; Woods et al.l 
1996). It has been shown that thermal winds arise at 
0.1-0.2 Rc, the Compton radius. Rc is defined to be 
R c ~ 10 10 (M/M Q )T C ; 8 1 cm, wher e T cs is the Comp - 
ton temperature in terms of 10 8 K (|Woods et al.i ri996). 
Therefore, to launch a thermally driven wind, we require 
the launching radius, -Riaunch, be located at greater than 
0.1 Rc- If we then assume that the observed velocity is 
equal to the local escape velocity, we can solve for the 
corresponding radius, 
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Setting this radius to be greater than or equal to the 
launching radius, which is required if the wind is to be 
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thermally driven, we find, 
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We can see here that for typical velocities, and low 
Compton temperature (< 10 6 K), driving winds by ther- 
mal pressure is difficult. We note that the actual veloc- 
ity is likely to be greater than the line-of-sight velocity 
due to inclination effects and transverse velocities across 
our line-of-sight. If the observed velocity is proportional 
to the gravitational potential, i.e., t> 2 ut > GM/i?i auncn , 
then a higher velocity would place the gas deeper in the 
potential well and thus increase the temperature needed 
to launch a thermally driven wind. In addition, we have 
assumed that the launching radius is the radius at which 
the observed velocity equals the escape velocity. 

Requiring that the velocity exceeds the escape velocity 
also requires that the bolometric luminosity to be 



Lnoi > (6A)-V\R/Rc)- 1/2 L C r 



(16) 



where Lqr is the critical luminosity defined as Lcr = 
2.88 x 10- 2 T~ 8 1/2 L E dd- (See lProgaetalll2000L for more 
details.) If the source luminosity is LboI < 2 x 10~ 2 Z/Edd 
then it would fail to launch an escaping wind. As seen 
in Figure [5j a majority of the wind sources are above 
this threshold, so thermal driving is plausible as long as 
our assumption about the launching radius is correct. If 
the wind originates closer than 0.1 Rq, then other mech- 
a nisms are needed. 

iLuketic et al.l (|2010t l perform hydrodynamical simula- 
tions to explore whether thermal driving could be respon- 
sible for the winds seen in X-ray binaries. They conclude 
that at low densities, thermal driving is possible from 
an X-ray heated accretion disk. However, at densities 
higher than n e > 10 12 cm" 3 , Compton heating is not 
sufficient at driving winds at velocities of v out > 10 2 km 



s . ILuketic et al.1 <|2010f ) compare their work to obser- 
vations of GRO J1655-40 which ha s a high density of 
n e ~ 10 14 cm" 3 (jMiller et al.ll2008t ). and conclude that 
thermal driving is not responsible for its winds. It is pos- 
sible that the other X-ray binaries have similar densities, 
and as Figure [1] shows, they have similar velocities as well 
as high ionizations. Therefore, Compton heating may be 
an unlikely driving source for these X-ray binaries. 

The AGN in Figure Q] show a much wider range in 
ionization but are all outflowing at velocities consis- 
tent with v ou t > 10 2 km s . If they share a simi- 
lar density to that of the BHB, then Compton heat- 
ing is not a viable drivin g mechanism for them either. 
iDorodnitsvn et all (|2008f ) also perform hydrodynamical 
simulations for Compton heated and radiation driven 
AGN winds. They find similar outflowing velocity and 
ionization parameters as we show in Figure[T]but the den- 
sity they assume is far lower than what is inferred from 
obser vations (e.g . . iMcKernan et al.ll2007t iRisaliti et all 
I2009t King et all |2012a|). In addition, the location of 
their warm-absorbers are much farther from the cen- 
tral source than t hose inferred from observat ions (e.g., 
IKing et al]|2012at ICrenshaw fc Kraemerll2012[ ). 



A third driving mechanism can be magnetic 
fields, whether through magneto-centrifugal force 
(Blandfor d fc Pavnel I1982D . or magnetic pressure from 
the toroid al field genera t ed by MRI in the disk as su. 
gest ed bv iContopoulosI (|1995[ ). iMiller fc Stonet (|20O 
and iProgal (|2003[ ). These winds tap magnetic field 
energy generated or sustained in the disk. It has 
been shown that for at least three of the sources in- 
cluded in our study, GRO 1655-40, NGC 4051, and 
NGC 4151, ma gnetic processes are likely driving the 
observed w i nds (iKraemer et al.l 120051: [Miller et all 120081 : 
IKing et al.l l2012at iNeilsen fe Homanl 120121 ). Because 
these sources (the BHB, GRO 1655-40, and Seyfert-l's, 
NGC 4051 and NGC 4151) are included in this relation, 
and span orders of magnitude in mass, it raises the possi- 
bility that magnetic forces may drive this wind relation. 

Jets are also thought to be driven by mag netic pro- 
cesse s in the disk or near the bla c k hole fe.g.. [Lovelace 
1971 IBlandford fc Pavnel [19821 : iBlandford fc Znaiek 
19771: iKrolik fc Hawlevl l20l6Th The jet power re- 
lation is determined to be log(Lj c t 42) = (1.18 ± 
0.24) log(L Bo i,42) - (0.96 ± 0.43). If we exclude Cygnus 
X-l, then the jet relation becomes log(Lj c t 42) = (1.34 ± 
0.50) log(i B oi,4 2 ) - (0.80±0.82) (See Figure©. When we 
include Cygnus X-l in the jet relation the slope of the jet 
relation and the initial wind relation are inconsistent at 
the 1.6cr level. When we exclude Cygnus X-l, the jet re- 
lation slope is formally consistent with that of the single 
wind relation. Further, if we examine the individual wind 
fits, the jet relation including or excluding Cygnus X-l is 
consistent with the BHB slope. When comparing the jet 
relation to the AGN sample, the slopes are only inconsis- 
tent at the 1.4cr and 1.2a level, including and excluding 
Cygnus X-l, respectively. Finally, we very tentatively 
suggest that the plane fit to the wind sample is consis- 
tent with the fundamental p lane of accretion onto black 
holes (|Gultekin et alj|2009 a). which would be further ev- 
idence of the similar dependence on mass accretion rate 
(as well as mass) of both the winds and jets. 

If the two types of outflows are regulated by the mass 
accretion rate in the same fashion, then the same driving 
mechanism may also be at work, as the geometry or mass 
loading of the magnetic fields may be driven by the mass 
accretion rate as well. This may explain the formal con- 
si stency of the slopes betwee n the jet and wind relations. 
A lBlandford fcTav nc (1982) scenario may be a viable so- 
lution fo r driving these outflows, an d could be possibly 
aided by Blandfor d fc Znaie scenario for jets. In 

addition. lOhsuga fc Mineshigd (|201lD show that MHD 
accretion flows can drive both jets and winds depending 
on the mass accretion rate, qualitatively consistent with 
Figure El 

The ultra-fast outflows appear to follow the jet rela- 
tion (See Figure [6|). These are winds whose observed 
velocity exceeds v > 0.01c. This raises the question of 
how these ultra-fast winds are accelerated to such high 
velocities. Are we seeing the phase at which these winds 
are being collimated into jets, as the power associated 
with the winds is very comparable to the jet power? 
Again, this could point to a shared driv ing mechanism 
between winds and jets, such as MHD (iLovelacd 119761 : 
IBlandford fc Pavnel fl982T) . if we are truly observing this 
transition phase between the two. 
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5.2. Implications for Feedback 

The characterization of these outflows allows us to de- 
termine that X-ray AGN winds are more efficient at re- 
movin g material than a re X-r ay BHB winds. Interest- 
ingly, Hopkins fc Elvisl (|2010l) show that only 0.5% of 
the bolometric luminosity needs to be converted into me- 
chanical power in order to regulate black hole growth and 
affect feedback in the host galaxy. As shown in Figure [TOl 
the majority of the AGN lie above (or are consistent 
with) 5 x 10~ 3 £bo1- A few sources lie above 5 x 10 ~ 2 I/R n i 
(the dotted line). iCrenshaw fc Kraemerl (|2012t ). in a 
study focusing only on AGN winds, show that up to half 
of their AGN are consistent with > 5 x 10~ 3 £bo1- How- 
ever, if the filling factor is much less than unity, the wind 
power will be far less than the 5 x 1CP 3 £bo1 limit for in- 
fluential feedback. This may imply that the X-ray winds 
do not have a large impact on feedback. 

In addition, Figure [5] shows a stratification of the ki- 
netic wind luminosity as a function of ionization in the 
AGN. The low ionization components (log£ < 2) tend to 
have a much higher kinetic luminosity as compared to the 
medium ionization components (2 < log£ < 3) and high 
ionization components (3 < log£). The reason for this 
may again be because the filling factor is not included in 
this analysis. As mentioned in Section|3l these black hole 
X-ray winds are thought to be clumpy and filamentary. 
Moreover, observations of ionized stellar winds indicate 
that the less ionized gas should be more clumpy, i.e. have 
a lower filling factor, due to pr essure confineme nt from 
the hot surrounding gas (e.e.. iSako et al.l[l999t ). This 
would imply that the low ionization components seen in 
Figure [5] would likely have a lower filling factor than the 
high ionization components. If the filling factors were 
included, the low ionization components would no longer 
rise above the higher ionization components. This would 
serve to flatten the initial wind relation, making the wind 
slope even more consistent with the jet relation slope. 
Feedback from the lowest ionization components would 
no longer dominate the relation. 

On the other hand, the low ionization components 
may also be consistent with being radiatively driven, 
and therefore would not follow the same relation as the 
high ionization components anyway. As shown in Section 
14. 1[ the high ionization components do follow a shallower 
slope of a Io s? >2 = 1.42 ± 0.06, which is consistent with 
both jet relations, i.e. including and excluding Cygnus 
X-l. ' 

Regardless of the ionization of the winds, jets are more 
efficient at a given bolometric luminosity, compared to 
X-ray winds. When considering power alone, jets may 
have a greater impact on mechanical feedback and galaxy 
evolution then winds. Depending on the mass accretion 
rate, for which we use the bolometric or Bondi luminos- 
ity as a proxy, we can now characterize the associated jet 
and wind power. Figure [5] shows exactly how the outflow 
power scales with Eddington fraction. There is a division 
between dominant outflow at approximately 10~ 2 LEdd- 
If both the winds and jets share a common launching 
mechanism, this division may be strongly driven by the 
mass accretion rate. Mass loading or even the geometry 
of the magnetic fields in the disk would have an impor- 
tant role as well, and can again be directly regulated by 
the mass accretion rate. The transition seen at approx- 



imately 10~ 2 LEdd is also interesting because this is the 
regime where winds begin to prevail over jet production, 
especially seen the spectral state dependence in X-ray 
binaries. 

We can now describe the outflow power as a function 
of Eddington fraction directly associated with the inner- 
accretion disk surrounding a black hole, vital for cos- 
mic simulations. This is important because as Figure [5] 
demonstrates, outflows are present in a range of Edding- 
ton fractions, not just low Eddington fractions. As galax- 
ies evolve through their "Quasar" and "radio" modes of 
accretion, we are still able to prescribe the outflowing 
power to assess the mechanical feedback in those systems 
and explore the implications for galactic co-evolution. 

5.3. Potential Caveats 

Before using these descriptions, it is important to 
understand the caveats involved in assembling this 
data set. As shown in Figure [SJ there is a poten- 
tial state dependence of out flow type on accretio n rate 
(also see IMiller et al.|[2006al l200l fNeilsen fc Leil200l 
iPonti et al.ll2012f ). However, it is not clear if this is a 
result of a selection bias toward high luminosity AGN. 
One could imagine that at low X-ray luminosity, i.e. el- 
lipticals and BHB in the "low/hard" state, detections 
of winds could be hampered by low signal-to-noise. This 
would be most pertinent to our jet sample, which is dom- 
inated by low luminosity AGN accreting at low accretion 
rates. However, even if winds were to coexist in these low 
accretion rates, j ust as they do in Seyfert 1 NGC 4051 
(|King et al.ll2011l) . jet power is likely to dominate by or- 
ders of magnitude, as the wind power is proportional to 
the ionizing luminosity, which would be small. In addi- 
tion, in BHB strong limits to wind detection have been 
made in the "low/hard" state (e.g.. iNeilsen fc Leell2009H 
as well as strong upp er limits to jet pr oduction in the 
"high/soft" state fe.g. JKing eFaIll2012bD . Therefore, we 
stress that this work is focused on the dominant outflow. 

We also note the difficulty in placing upper limits on 
wind detections using absorption features in the X-ray 
band. As these features can be seen at an array of dif- 
ferent velocities and ionization states, there is no specific 
wavelength one would expect to find an absorption fea- 
ture denoting an outflow. Further, the wind power es- 
timates do not depend on the strength of the line, but 
only the wavelength and ionization state, making esti- 
mates of upper limits rather difficult. These issues of 
detection affect both BHB and AGN in the same man- 
ner, and we stress that the lack of detection of these 
absorption features is not evidence for the absence of 
a wind, but may be the absence of evidence. Again, 
the state dependence of outflows seen in BHB and now 
in AGN (Figure [5]) is likely to be driven by accretion 
rate. However, longer integrations to improve signal- 
to-noise of BHB in the "low/hard" state and AGN at 
low accretion r ates are needed to be confident of this as- 
sessment (e.g., IMiller et a.1.1 l2012h . Next, we note that 
there are outflows other than the ones examined in our 
analysis that are still important in removing substan- 
tial amounts of material from their accretion disks and 
host galaxies. In particular, broad absorption line (BAL) 
quasars have particular l y powerful outflows ()Moe et al.l 
[20091: iDunn et al.l l2010t iBrandt et ail [2000h . However, 
these outflows are observed in the optical and ultra-violet 
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regime and have much lower ionization parameters than 
the X-ray winds discussed here. Therefore, they are not 
as readily associated with - or driven by - the inner ac- 
cretion disk, and have not be en included in our analysis. 
ICrenshaw fe Kraemerl (|2012l ) also show a positive corre- 
lation between ionization parameter, U, and column den- 
sity in local AGN in their figure 3. This demonstrates 
that the bulk of the outflow material is being observed 
in the X-ray regime. 

A broader range of ionization parameters are probed in 
the AGN as compared to the BHB (See Figure [TJ, which 
may also contribute to the AGN scatter. The mix of £ in 
AGN calls into question whether we are probing the same 
physics, i.e. closest to the black holes. Because ioniza- 
tion is dependent on the distance as £ = L(nr 2 ) -1 , simi- 
lar ionization states should probe the same distance from 
the black hole for a given luminosity. Consequently, sim- 
ilar micro-physics at a given radius and ionization should 
be at work. Further, in Section 15.21 examining the high 
ionization components alone results in a shallower slope 
when comparing bolometric luminosity to wind power. A 
shallower relation is more consistent with the jet relation. 
This demonstrates the clear need for a much larger sam- 
ple size. Fortunately, Astro- H will provide the needed 
coverage in the highest ionization band. This will not 
only allow for the detection of additional sources, but 
also detection of the highest ionization states for com- 
parison with stellar-mass black holes. 

Although we see a large range in ionization, we do 
not see as large of a range in velocity. It is only when 
we include the ultra-fast outflows that three orders of 
magnitude in velocity are probed as compared to the six 
orders of magnitude in the ionization parameter. This 
trend is important in understanding whether inclination 
has an effect on the given L W md vs LboI correlation. As 
BHB winds are thought to be o bs erved in nearly ed ge-on 
sources (e.g. IMiller et al.l l2006a b; P onti et al.ll2012T) . and 
AGN winds, especially Seyfert 1 AGN winds, are though t 
to be observed in face-on sources (e.g. JWu fe Haai200l . 
inclination could have the potential to bias our results. 
However, the data show no trend in velocity as a func- 
tion of inclination. In addition, when examining BHB 
sources individually, face-on sources (e.g., GX 339-4, 
XTE J1817-330) do not show absorption features in the 
Fe K band when they are in the "high/soft" state, con- 
trary to their edge-on counterparts. This is due to lim- 
ited sensitivity, since low inclination sources tend to be 
softer a nd to give less signa l through the Fe K band. As 
noted in iPonti et al.l (|2012f) , the limits on flux in face-on 
sources are not very constraining, and lines as weak as 
those in H1743-322 (|Miller et al.ll2006bft could not have 
detected in e.g. XTE J1817-330. For instance, there 
is likely a simple absence of evidence for BHB winds in 
face-on systems. 

One may also expect inclination to have an effect on the 
estimated kinetic jet power. For example, those sources 
for which the jet is directed along our line of sight may 
suffer from Doppler boosting. However, this would pri- 
marily affect the radio luminosity of such sources, and 
not the kinetic power, which is taken from estimates of 
cavity sizes. On the other hand, the jet power may be 
influe nced by the spin of the black hole. If jets are driven 
by the iBla ndford fe Znaiekl (|1977l ) mechanism, then spin 
may play a large role in the power released by the jets. 



However, the common slope s between the jet and wind 
relation points to more of a IBlandford fe Pavnd (fl982^ 
scenario, where the spin of the black hole does not af- 
fect the power released. Further, the fact that the jet 
and wind power seem to be present at certain Eddington 
fractions, point to the idea that mass accretion rate may 
be the throttle that is ultimately driving the type and 
power of the outflow. We note that spin is unlikely to 
play a large role in the X-ray wind power regardless, as 
winds are thought to originate further out in the accre- 
tion disk. 

An additional concern with the jet power is that the 
estimate is a long-term average and not instantaneous 
as are the wind power estimates. If the jet production 
occurs on timescales that are much shorter than the dy- 
namical timescale of the cavity, than the power estimates 
would increase. Unfortunately, this is a limit of this tech- 
nique when using cavities to estimate power. However, 
long term estimates of power are more pertinent for feed- 
back estimates. 

Finally, as this sample is small in size, the results must 
be regarded cautiously and tested in the future. It is im- 
perative that we obtain more observations at all masses 
and mass accretion rates. Specifically, black holes ac- 
creting at LboI ~ 10 41 — 10 42 ergs s _1 could distinguish 
whether one linear fit is required across the entire wind 
sample or if the BHB and AGN are better fit by individ- 
ual linear fits. This could either be a small Seyfert galaxy 
with mass on order of M ~ 1O 5 M accreting at a few 
percent of Eddington, or a large SMBH, M ~ 10 9 M©, 
accreting at a very low Eddington rate. In addition, 
non-simultaneity of AGN luminosities could have a dra- 
matic effect on the observed scatter seen in the X-ray 
winds. Although the AGN timescales for disk evolution 
are longer than BHB, observations made years apart may 
not probe the same accretion regime. 

6. CONCLUSIONS 

• In this study, we find that winds are consistent 
with being regulated according to a simple re- 
lation across a large mass scale. In particular, 
we find the trend is described as log L W ind,42 °c 
(1.58 ± 0.07) logLBoi,42- The slope is greater than 
unity, so it may imply that the SMBH are more 
efficient at expelling material than BHBs. 

• If we fit the BHB and AGN populations separately, 
they still require consistent slopes of a BHB = 
0.91 ± 0.31 and a AGN = 0.63 ± 0.30. Further, if 
we assume mass is influencing this relation and fit 
a plane to the data, we find the best fit relation 
to be log(L wind ) = (1.2 ± 0.3)log(M BH ) + (0.2 ± 
0.4) log(L Bo i) + (24.5 ± 0.2) with scatter cr = 0.68 
consistent with the "fundamental plane" of accre- 
tion onto black holes. 

• It remains possible that different processes tied to 
the mass accretion rate- thermal driving in stellar- 
mass black holes and radiative driving in AGN - 
are actually at work in driving winds. However, it 
is not clear that these different mechanisms should 
agree so well and follow the same slope in these 
wind relations. Moreover, it seems that a mag- 
netic wind must be at work in GRO 1655—40, NGC 
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4051 a nd NGC 4151 (iMiller et al.l200 8 HKing et all 
I2012at iKraemer et al l 120051 ). which fall on the re- 
lation. This may also suggest a role for magnetic 
driving across the mass scale. 



• Furthermore, when we examine jet power, the 
data may be consistent with winds and jets be- 
ing regulated in a common fashion. Since radia- 
tive and thermal processes are not likely to drive 
relativistic jets, a mechanism like magnetocentrifu- 
gal or MHD winds are plausible expla nations (e.g., 
IBlandford fc PavndflQSl IProgall200lh . 

• The ultra-fast winds appear to obey the same reg- 
ulation scheme as slower, more common winds, if 
they have a low filling factor and the slow winds 
have a high filling factor close to unity. However, 
some ultra-fast winds appear to carry as much ki- 
netic luminosity as jets, even after accounting for 
filling factors. This suggests that we may be see- 
ing a phase where winds finally are accelerated into 
jets. 



• Figure [5] provides a direct way to quantify the out- 
flow power as a function of mass accretion rate. 
A division between dominant outflow state is ob- 
served at approximately 10 _2 LEdd- This trend has 
broad implications, especially for theoretical sim- 
ulations that need prescriptions for feedback to 
study galactic dynamics and evolution. 

• A larger sample will help us to distinguish between 
these proposed relations as well as quantify the in- 
trinsic scatter. As it stands now, Chandra will play 
an integral part in future studies. Looking further 
ahead, Astro-H will have improved sensitivity in 
the Fe K band, enabling unprecedented looks at 
the most ionized and innermost flows in the accre- 
tion disks of both BHB and AGN. 
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Individual X-ray Wind Components 
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NGC 3516 


SI 


1 


2.40 ± 0.15 


r i r t'\ 1 "1 A T 

950 ± 147 


41.55 ± 0.26 


McKcrnan ct al. (2007) 


JNGC O/OO 


C 1 

ol 


1 


o nn _L n i n 
2.yo ± U.1U 


£n £ _l_ i £ 
505 zt 15 


a 1 no _L n i o 
41.02 it 0.12 


^ I,.. I," i. ...I ion^7l 

McKcrnan ct al. ( zOU / ) 






2 


2.10 ± 0.10 


515 ± 15 


41.85 ± 0.12 






3 


0.40 ± 0.10 


545 ± 25 


43.62 ± 0.12 








4 


3.00 ± 0.10 


1145 ± 42 


41.99 ± 0.12 




NGC 4051 


SI 


1 


1.00 ± 0.30 


520 ± 82 


41.26 ± 0.37 


McKcrnan et al. (2007) 






2 


2.60 ± 0.25 


600 ± 77 


39.85 ± 0.30 








3 


3.80 ± 0.10 


2230 ± 55 


40.36 ± 0.11 




NGC 4051 


SI 


1 


4.50 ± 0.90 


680 ± 40 


37.81 ± 0.91 


Kine et al. (2012b) 






2 


3.28 ± 0.04 


640 ± 45 


38.95 ± 0.11 






3 


1.00 ± 0.11 


400 ± 325 


40.62 ± 1.06 




NGC 4151 


SI 


1 


3.58 ± 0.30 


491 ± 8 


40.30 ± 0.30 


Kracmer ct al. (2005) 






2 


2.26 ± 0.30 


491 ± 8 


41.62 ± 0.30 


NGC 4593 


SI 


1 


2.61 ± 0.90 


400 ± 121 


40.47 ± 0.98 


McKcrnan et al. (2007) 






2 


0.50 ± 0.30 


380 ± 137 


42.51 ± 0.56 


NGC 5548 


SI 


1 


2.20 ± 0.20 


560 ± 77 


41.77 ± 0.27 


McKcrnan et al. (2007) 






2 


3.90 ± 0.15 


830 ± 172 


40.59 ± 0.31 


BHB 


4U 1630 




1 


4.90 ± 0.40 


300 ± 200 


32.68 ± 0.96 


this paper 


GRO 1655-40 




1 


4.90 ± 0.20 


500 ± 200 


32.31 ± 0.56 


Miller ct al. (2008) 


GRO 1655-40 




1 


4.20 ± 0.15 


470 ± 230 


33.42 ± 0.66 


Ncilscn & Homan (2012) 


H 1743-322 a 




1 


5.50 ± 0.10 


670 ± 170 


33.43 ± 0.35 


this paper 


H 1743-322 b 




1 


5.60 ± 0.10 


340 ±170 


32.32 ± 0.66 


this paper 


GRS 1915+105 si 




1 


4.30 ± 0.20 


300 ± 200 


33.42 ± 0.89 


Ucda ct al. (2009) 


GRS 1915+105 s2 




1 


5.60 ± 0.20 


1000 ± 200 


34.09 ± 0.33 


Miller ct al. (2006b) 


GRS 1915+105 s3 




1 


5.50 ± 0.50 


1400 ± 300 


34.35 ± 0.57 


this paper 


GRS 1915+105 s4 




1 


6.00 ± 0.40 


1100 ± 400 


33.68 ± 0.62 


this paper 



GUIS 1915+105 s5 I 6.20 ± 0.70 900 ± 400 33.22 EE 0.91 this paper 



Note. — The above table lists all the components that are considered in this analysis. SI stands for Scyfcrt 1 and S2 
stands for Seyfert 2. The kinetic luminosity from the AGN that have more than one component for a single observation 
are summed and included as total kinetic luminosities in Table [2] 
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TABLE 2 
X-ray Wind and Jet Quantities 







Object 


log L B oi 


log L wint i/Cv 


log M B h 


log£> 


Code 


Reference 








(ergs 




(ergs s — 1 ) 


(M e ) 


(cm) 






X-ray winds 




1 


Akn 564* 


44.50 


± 0.13 


42.52 ± 0.63 


6.9 e 


26.51 


XSTAR 


McKcrnan ct al. (2007) 




2 


IC 4329a* 


43.80 


± 0.13 


42.15 ± 0.85 


7.0 C 


26.32 


XSTAR 


McKernan et al. (20071 




3 


IRAS 18325* 


44.60 


± 0.13 


41.98 ± 0.36 


7.0 6 


26.42 


XSTAR 


Zhang et al. (2011a') 




4 


NGC 3516 


43.50 


± 0.13 


41.55 ± 0.26 


7.5" 


26.07 


XSTAR 


McKernan et al. (2007') 




5 


NGC 3783 


44.20 


± 0.13 


43.64 ± 0.12 


7.5 C 


26.11 


XSTAR 


McKernan et al. (2007) 


AGN 


6 


NGC 4051 


42.60 


± 0.13 


41.32 ± 0.32 


6.2 a 


25.50 


XSTAR 


McKernan et al. (2007) 


7 


NGC 4051 


42.60 


± 0.13 


40.63 ± 1.04 


6.2 a 


25.50 


XSTAR 


King et al. (2012aT" 




8 


NGC 4151* 


43.90 


± 0.13 


41.64 ± 0.29 


7.1 c 


25.64 


XSTAR 


Kraemer et al. (20051 




9 




43.70 


+ 0.13 


49 _|_ fl cr 
HI KJ.OO 


6.7 C 


26.04 


XSTAR 


McKorTian pt al (9007^ 




10 


NGC 5548 


44.30 


± 0.13 


41.80 ± 0.26 


7.6" 


26.35 


XSTAR 


McKernan et al. (2007) 




11 


MCG -6-30-15* 


43.40 


± 0.13 


40.89 ± 0.22 


6.5 d 


26.01 


XSTAR 


McKernan et al. (20071 




12 


Mrk 290 


44.40 


± 0.13 


42.57 ± 0.36 


7.4 a 


26.60 


Cloudy 


Zhang~al. (2011b^ 




13 


Mkn 509 


45.20 


± 0.13 


41.93 ± 0.54 


8.2 C 


26.65 


Cloudy 


Ebrcro ct al. (2011) 




14 


4U 1630 


38.20 


± 0.43 


32.68 ± 0.96 


1.0* 


22.42+0.30-? 


XSTAR 


this paper 




15 


GRO 1655-40a 


37.70 


± 0.29 


32.31 ± 0.56 


0.83 9 


21.79+0.20' 


XSTAR 


Miller et al. (20081 


BHB 


16 


GRO 1655-40b 


37.80 


± 0.29 


33.42 ± 0.66 


0.83 9 


21.79±0.20 i 


XSTAR 


Neilsen & Homan (2012) 


17 


H 1743-322 a 


38.60 


± 0.46 


33.43 ± 0.35 


1.0* 


22.42+0.30-'' 


XSTAR 


this paper 




18 


H 1743-322 b 


38.50 


± 0.43 


32.32 ± 0.66 


1.0'' 


22.42+0.30-? 


XSTAR 


this paper 




19 


GRS 1915+105 SI 


38.90 


± 0.13 


33.42 ± 0.89 


1.15 h 


22.54±0.03 m 


XSTAR 


Ueda et al. (20091 




20 


GRS 1915+105 S2 


39.50 


± 0.13 


34.09 ± 0.33 


1.15 h 


22.54±0.03 m 


XSTAR 


Miller et al. (2006b1 




21 


GRS 1915+105 S3 


39.00 


± 0.13 


34.35 ± 0.57 


1.15 h 


22.54±0.03 m 


XSTAR 


this paper 




22 


GRS 1915+105 S4 


39.10 


± 0.13 


33.68 ± 0.62 


I.15 h 


22.54±0.03 m 


XSTAR 


this paper 




23 


GRS 1915+105 S5 


39.10 


± 0.13 


33.22 ± 0.91 


I.15 h 


22.54±0.03 m 


XSTAR 


this paper 


jets 






l°g J^Bondi 


log L Jet 
















(ergs 


s-i) 


(ergs s" 1 ) 












24 


NGC 507 


44.41 


± 0.09 


44.01 ± 0.15 


8.9' 


26.34 




Allen et al. (20061 




25 


NGC 4374 


43.69 


± 0.30 


43.18 ± 0.13 


8.8 f 


25.72 




Allen et al. (20061 




26 


NGC 4472 


43.79 


± 0.25 


42.91 ± 0.13 


8.9f 


25.72 




Allen et al. (20061 




27 


NGC 4486 


44.16 


± 0.35 


43.54 ± 0.23 


9.5 f 


25.72 




Allen et al. (20061 




28 


NGC 4552 


43.37 


± 0.22 


42.19 ± 0.11 


8.7 f 


25.72 




Allen et al. (20061 


AGN 


29 


NGC 4636 


42.29 


± 0.24 


41.48 ± 0.12 


8.2? 


25.72 




Allen et al. (20061 




.30 


NGC 4696 


43.40 


± 0.56 


42.90 ± 0.17 


8.6 f 


26.14 




Allen et al. (20061 




31 


NGC 5846 


42.85 


± 0.42 


41.87 ± 0.16 


8.6? 


25.88 




Allen et al. (20061 




32 


NGC 6166 


43.49 


± 0.30 


43.20 ± 0.13 


8.9 f 


26.62 




Allen et al. (20061 


BHB 


33 


Cygnus X-l 


37.30 


± 0.13 


36.00 ± 2.00 


1.0' 


21.76±0.04 fe 




this paper 


Ultra-Fast Outflows 




34 


3C 111 


45.9 




45.8 


9.0 


26.81 


XSTAR 


Tombcsi ct al. (2011) 


AGN 


35 


APM 08279+5255 


47.3 




50.5 


10.3 


28.3 




Chartas et al. (20021 




36 


PC 1115+080 


46.5 




51.2 


9.0 


28.1 




Chartas ct al. (2007) 


BHB 


37 


J17091+3624 


37.5 




38.3 


1.0* 


22.41 


XSTAR 


King ct al. (2012b) 



Note. — * These sources have bole-metric luminosities estimated from their 2-10 kcV fluxes ( L Boi 20L 2 -iofc e v , Vasudcvan & Fabian 2009 ). Th e 
m asses are give n by e ach reference unless otherwise stated; a : iDenney et ah] 120HJ), b: Lcc (2005), c: [Peterson et alj (20041) , d : IMcHardy et all 420051) , 
e: IColl ier et all i2001l) . /: derived from the relation given in Trcmainc ct al. (2002) using the a given in lAllen et al 1 ll2006D . g: IShahbaz et al.l (Il99gt )7fa 
IGreiner et al.1 i2001l) . i: the mass of these BHB has no t been det ermined so a M — IO Mq and 20% error has bee n assumed j: the distance to these sources is 
unknown and assumed to be 8.5+4 kpc, k: IReid et all 120111) . /: IFoellmil 02009) , m: IHarlaftis fc Greinerl 1 120041) 
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TABLE 3 
Individual X-ray Wind Components 



Data Set a ft 7 ctq 



Winds 



ALL 


1.58±0.07 


-(3.19±0.19) 


0.68 


BHB 


0.91±0.31 


-(5.58±1.68) 





AGN 


0.63±0.30 


-(1.24±0.63) 


0.58 


log? > 2 


1.42±0.06 


-(3.73±0.14) 


0.56 


ALL 


0.2±0.4 


1.2±0.3 


24.5±0.2 0.68 


Jets 


ALL 


1.18±0.24 


-(0.96±0.43) 





AGN 


1.34±0.50 


-(0.80±0.82) 






Note. — These are the best fit parameters for each of our 
linear models. The a. parameter describes the coefficient of 
the bolomctric or Bondi luminosity, j3 is the normalization 
of each linear fit except for the last wind fit. In that case it 
is the coefficient of the mass term and 7 is the normalization 
of the fit. Finally, ctq is the intrinsic scatter of each fit. 



s BHB 

1 5 T AGN 

: ufo 



W 1Q 4L 



E 



8 10 3 

> 



10' 



i — Hf 1 



- , 1 , , L 



4 6 

log (0 



Fig. 1. — The above figure plots the observed velocity components versus ionization, for the slow and "ultra-fast" winds in our black hole 
sample. In black are the AGN winds, in blue are the BHB winds and in red are the ultra-fast winds. The points with arrows denote lower 
limits to the ionization state, as the actual state for these ultra fast winds was not analyzed with a photoionization model. 
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-4-2 2 

lo g( L boi,4 2 ) 



Fig. 2. — The plot above shows the correlation between bolometric luminosity and kinetic wind luminosity in individual outflowing 
components. The line represents the best fit to the total kinetic luminosities which are plotted in Figure [3] while the yellow line is the 
best fit to the individual components with log£ > 2. The high ionization parameters are described by the following form log(L w ; nt i 42) = 

(1.42 ± 0.06) log(L Bo i,42) ~ (3.73 ± 0.14), with an intrinsic scatter of a l g g6>2 = 0.57 



O 



O 




-2 

'ogC-boi^a) 



Fig. 3. — The plot above shows the correlation between bolometric luminosity and kinetic wind luminosity. The black line is described by 
log(L w ; n( j 42) = (1.58±0.07) log(Z/Bol,42) — (3.19±0.19), with an intrinsic scatter of oq = 0.68. The blue dashed region is the lcr confidence 
region including the scatter of the relation. The solid region is the la confidence region excluding the scatter. The wind kinetic luminosity 
is plotted per filling factor. The plot shows a simple regulation of wind production across a large mass scale, and the slope indicates that 
the SMBH winds are more efficient then the stellar-mass black holes. 
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lo g( L boi,4 2 ) 

Fig. 4. — The plot above depicts the wind power versus the bolometric luminosity, just like Figure \3\ but this figure includes the jet 
power as red data points. The red line describes all the jet points as log(Lj ct 42) = (1.18 ± 0.24) log(L Bondi 42) — (0.96 ± 0.43). The yellow 
line describes the data set if Cygnus X-l is excluded from the fit is given as log(Lj et 42) = (1.34 ± 0.50) log(LBondi,42) — (0.80 ± 0.82). 
The dashed regions are the la confidence regions. The orange line and dashed region is the best fit line and la confidence region when 
excluding Cygnus X-l from the fit. One can see that the normalization of the jets is higher, demonstrating that for a given bolometric 
luminosity they are more powerful. One can also see that the slope between the two relations is quite similar, perhaps indicating a common 
launching mechanism 




Fig. 5. — The above plot shows the power emitted either from the jet power (AGN:red and BHB:orange) or wind power (AGN:black and 
BHB:blue) as a compared to the mass accretion rate, which is approximated by the bolometric luminosity on the x axis. A clear turnover 

at Mace ~ 10 _2 MEdd indicates where the power emitted is becoming less efficient. Interesting is the dichotomy between where the jets lie 
at lower mass accretion rates and where the winds lie at higher accretion r ates. The thick black li ne denotes the output power by outflows, 
where as the thin line is the power generated by radiation as described by Churazov ct al. (2005) 



22 A. L. King et al. 




-4 -2 2 4 6 
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Fig. 6. — This plot is the same as Figure[4]but now includes UFO's in black squares (v > 0.1c). The upper square is the power estimate 
with a filling factor of unity. The bottom square connected by the dashed line is the lower estimate of the wind power if the filling factor is 
as low as C v = 10 -4 . Even with a smaller filling factor, the UFO's resemble the jet relation more so then the wind relation. Perhaps this 
is indicated that the winds are reaching a phase where they are being accelerated into jets. 




-4 -2 2 4 

| og(L b oi,4 2 ) 



Fig. 7. — This figure shows the best fit linear regressions when the AGN and BHB samples are fit separately. The BHB are described 
by log(i^™ 42) = (°- 91 ± °- 31 ) 1 °g( i IS B 42) ~ ( 5 ' 58 ± L68 ) with scatter consistent with zero. The AGN are described by log(L^f 42 ) = 
(0.63 ± 0.30) log(L^N 2 ) _ ( L2 4 ± 0.63) with scatter afi- GN = 0.58. 
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30 I i i i i i i i i i i i i i i i i i i i i i i i 

8 10 12 14 16 18 20 
1.22log(M BH ) + 0.22log(L Bol ) 

Fig. 8. — The above plot shows the best fit plane of our wind sample when including mass as a third parameter with bolometric luminosity 
and wind power. The plane is described by log(L willd ) = (1.2 ± 0.3) log(M B H) + (0.2 ± 0.4) log(L Bol ) + (24.5 ± 0.2) with scatter cr = 0.68. 




0.78log(M BH ) + 0.67log(L Bol ) 

Fig. 9. — This plot shows our wind data plotted against the fundamental plane of black hole activity described by Giiltckin et al. (2009a). 
The solid line is used to show the one-to-one correspondence in this plane cross section. t/=5 GHz. Although the inte rcepts are different, the 
coefficients of mass and X-ray/Bolometric luminosity are consistent between our sample and Giiltckin ct al. (2009a), which may tentatively 
suggest a common driving mechanism. 
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Fig. 10. — This plot shows the correlation between kinetic wind luminosity per filling factor divided by the bolometric luminosity as 
compared to the bolom etric luminosity. The do tted line is 5% I/ Rnl , while th e dashed line is 0.5% £bo1- These are the limits of the kinetic 
wind power reported by Di Mattco ct al. (2005) and Hopkins & Elvis (2010), respectively, for mechanical feedback to have an influence on 
black hole growth and feedback. We expect these winds to have a small filling factor which would make the wind power estimate lower, 
and perhaps below 0.5% LboI- 



